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PREFACE 

Preparation  of  the  NBS  Self-Study  Manual  on  Optical  Radiation  Measurements  is  an 
important  part  of  our  effort  to  meet  the  needs  of  the  National  Measurement  System  with 
respect  to  the  measurement  of  optical  electromagnetic  radiation.   Significant  needs  exist 
for  measurements  with  uncertainties  of  about  one  percent,  but  the  uncertainties  actually 
achieved  are  often  5   to   10  percent,   or  even  more.   These  relatively  poor  accuracies, 
as  compared  to  those  in  many  other  types  of  physical  measurements,  result  to  a  great 
extent  from  the  multidimensional  character  of  optical  radiation;  radiant  energy  is  dis- 
tributed and  may  vary  with  position,  direction,  wavelength,  polarization,  and  time. 
Measurement  results  are  also  affected  by  various  instrumental  and  environmental  parameters. 
In  addition,  many  of  those  who  make  measurements  of  optical  radiation  have  little  or  no 
training  or  experience  in  this  field  and  are  limited  in  the  amount  of  time  they  can  devote 
to  acquiring  the  needed  information,  understanding,  and  proficiency  with  measurement 
techniques.   Moreover,  there  are  few  schools  that  offer  courses  in  radiometry,  and  there 
are  almost  no  adequate  texts  or  references  dealing  with  this  entire  subject. 

The  idea  of  producing  a  Self-Study  Manual  at  NBS  to  try  to  fill  some  of  this  void 
was  developed  by  one  of  us  (HJK)  in  the  latter  part  of  1973.   Detailed  planning  got  under 
way  in  the  summer  of  1974  when  a  full-time  editor  (FEN)  was  appointed.   The  first  nine 
chapters  were  published  in  a  series  of  four  NBS  Technical  Notes  (910-  series)  up  to 
June  1979.   However,  subsequent  progress  was  slowed  by  administrative  uncertainties  and 
problems  that  affected  the  availability  of  authors  (causing  several  changes  in  scheduling) 
and  that  finally  ended  all  direct  NBS  funding  for  the  project  on  September  30,  1979. 
Welcome  transfers  of  funds  from  a  number  of  other  federal  agencies  were  obtained  to  carry 
the  project  to  September  30,  1980,  with  the  hope  that  it  could  then  be  restored  to  the  NBS 
budget.   When  that  failed  to  materialize,  outside  funding  was  again  obtained  to  Septem- 
ber 30,  1981  (through  FY  81). 

Meanwhile,  the  program  director  (HJK)  retired  at  the  end  of  May  1981,  and  the  editor 
(FEN)  on  September  29,  1981.   However,  remaining  FY  81  funds  were  used  for  a  grant  to 
Catholic  University  of  America  (CUA) ,  with  whom  NBS  negotiated  a  cooperative  agreement 
under  which  CUA  hired  the  full-time  editor  (FEN)  who  continues,  as  before,  at  NBS  (same 
address  and  telephone),  working  closely  with  the  new  project  head  (JBS)  there.   CUA  will 
also  engage  the  services  of  three  other  former  NBS  employees:   Dr.  Kostkowski  as  part-time 
consultant  throughout  the  project  and  as  author  of  two  more  chapters,  and  Dr.  Venable  and 
Mr.  Richmond  as  authors  of  one  chapter  each,  as  outlined  in  more  detail  below. 

This  Technical  Note  910-5,  with  the  first  chapter  of  Part  III — Applications,  is  a 
departure  from  Part  I — Concepts,  and  the  only  one  now  scheduled.   Its  objectives  and 
significance  are  explained  in  the  abstract  and  in  the  opening  section,  and  are  further 
amplified  below  in  a  brief  discussion  of  the  reasons  for  our  preoccupation  with  Part  I — 
Concepts. 

Our  aim  for  the  entire  Manual  has  been  to  provide  a  comprehensive  tutorial  treatment 
that  is  complete  enough  for  self  instruction.   That  is  what  is  meant  in  the  title  by 
"self-study";  the  Manual  does  not  contain  explicitly  programmed  learning  steps  as  that 
term  often  connotes.   In  addition,  through  detailed,  yet  concise,  chapter  summaries,  the 
Manual  is  designed  to  serve  also  as  a  convenient  reference  source.   Those  already  familiar 
with  a  topic  should  turn  immediately  to  the  summary  at  the  end  of  the  appropriate  chapter. 
They  can  determine  from  that  summary  what,  if  any,  of  the  body  of  the  chapter  they  want  to 
read  for  more  details. 

The  material  in  the  Manual  is  presented  at  the  level  of  a  college  graduate  in  science 
or  engineering,  but  even  for  those  with  facility  in  college  mathematics  and  a  first  course 
in  physics,  it's  not  at  all  easy  reading  in  spite  of  our  best  efforts  at  clarity  and  sim- 
plicity.  This  is  an  unavoidable  result  of  the  primary  aim  ("to  make  one-percent  measure- 
ments commonplace")  coupled  with  the  fact  that  it  must  serve  the  needs  of  so  many  different 
fields,  including  astronomy  and  astrophysics,  mechanical  heat-transfer  engineering,  illumi- 
nation engineering,  photometry,  meteorology,  photo-biology  and  photo-chemistry,  optical 
pyrometry,  remote  sensing,  military  infrared  applications,  etc. 
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Apparently  it  is  very  difficult  for  those  who  have  not  been  directly  involved  to 
realize  the  full  implications  of  that  situation.   Each  of  us  tends  to  think  of  radiometry 
and  radiometric  measurements  in  terms  of  our  immediate  experiences  and  requirements.   What 
each  of  us  would  like  to  have  is  a  set  of  simple,  carefully  designed  procedures,  with 
appropriate  cautions  concerning  likely  sources  of  error,  for  making  our  own  particular 
measurements.   However,  the  next  reader  wants  the  same  thing,  but  for  entirely  different 
measurements.   The  desired  radiometric  quantities  to  be  measured  are  different,  the  instru- 
mentation is  different,  the  ambient  conditions  are  different  —  the  possible  ways  in  which 
significant  differences  may  exist,  in  terms  of  the  radiation  parameters  (position,  direc- 
tion, spectrum,  time  or  frequency  of  modulation  or  fluctuation,  and  polarization)  as  well 
as  instrumental  and  environmental  parameters,  are  so  numerous  that  any  attempt  to  cover 
them  all  with  a  "cookbook"  treatment  of  specific  measurement  procedures  would  be  impossibly 
unwieldy  and  could  never  be  completed  within  any  conceivable  budget  limitations  short  of 
Utopia.   The  only  way  in  which  we  can  hope  to  effectively  assist  every  reader  who  needs  to 
make  one-percent  measurements  is  to  provide  you,  the  reader,  with  material  which,  with 
sufficient  effort  on  your  part,  will  help  you  to  develop  sufficient  understanding  and 
grasp  of  basic  principles  to  solve  your  own  particular  measurement  problems.   That's  why 
we  have  concentrated  on  the  basic  material  of  Part  I — Concepts  and  why  it  is  not  easy 
reading.   The  problems  of  radiometry  in  general  just  aren't  simple  ones,  even  when  limited 
to  classical  radiometry  in  the  domain  of  geometrical  optics  (ray  optics) .   Although  we 
have  kept  within  that  limitation  so  far,  we  recognize  the  increasingly  urgent  need  for 
similar  treatment  of  the  problems  of  coherent  radiation  and  we  are  happy  to  announce  that 
the  next  chapter  to  be  published,  in  a  few  months  we  hope,  will  be  Part  I,  Chapter  10, 
"Introduction  to  Coherence  in  Radiometry"  by  John  B.  Shumaker. 

The  basic  approach  and  focal  point  of  the  treatment  in  this  Manual  is  the  measurement 
equation,  first  introduced  in  detail  but  limited  to  the  radiation  parameters  of  position, 
direction,  and  spectrum,  in  Chapter  5.   We  believe  that  every  measurement  problem  should 
be  addressed  with  such  an  equation,  relating  the  quantity  desired  to  the  data  obtained, 
through  a  detailed  characterization  of  the  instruments  used  and  the  radiation  field 
observed,  in  terms  of  all  of  the  relevant  parameters.   These  parameters  always  include  the 
radiation  parameters  (listed  above) ,  as  well  as  environmental  and  instrumental  parameters 
peculiar  to  each  measurement  configuration.   The  objective  of  the  Manual  is  to  develop  the 
basic  concepts  required  so  that  the  reader  will  be  able  to  use  this  measurement-equation 
approach.   It  is  our  belief  that  this  is  the  only  way  that  uncertainties  in  the  measurement 
of  optical  radiation  can  generally  be  limited  to  one,  or  at  most  a  few,  percent. 

The  original,  overall  plan  for  the  Manual  organized  it  into  three  Parts: 

Part  I — Concepts: 

Step-by-step  buildup  of  the  measurement  equation   in  terms  of  the  radiation  parameters, 
the  properties  and  characteristics  of  sources,  optical  paths,  and  receivers,  and  the 
environmental  and  instrumental  parameters.   Useful  quantities  are  defined  and  discussed, 
and  their  relevance  to  various  applications  in  many  different  fields  (photometry,  heat- 
transfer  engineering,  astronomy,  photo-biology,  etc.)  is  indicated.   However,  discussions 
of  actual  devices  and  measurement  situations  in  this  part  are  mainly  for  purposes  of 
illustrating  concepts  and  basic  principles. 

Part  II — Instrumentation: 

Descriptions,  properties,  and  other  pertinent  data  concerning  typical  instruments, 
devices,  and  components  involved  in  common  measurement  situations.   Included  is  material 
dealing  with  sources,  detectors,  filters,  atmospheric  paths,  choppers  (and  other  types  of 
optical  modulators),  prisms,  gratings,  polarizers,  radiometers,  photometers,  spectroradiom- 
eters,  spectrophotometers,  etc. 

Part  III — Applications: 

Measurement  techniques  for  achieving  a  desired  level  of,  or  improving,  the  accuracy 
of  a  measurement.   Included  will  be  a  very  wide  variety  of  examples  of  environmental  and 
instrumental  parameters  with  discussions  of  their  effects  and  how  to  deal  with  them.   The 
examples  will  also  be  drawn  from  the  widest  possible  applications  in  illumination  engi- 
neering, radiative  heat  transfer,  military  infrared  devices,  remote  sensing,  meteorology, 
astronomy  and  astrophysics,  photo-chemistry,  photo-biology,  etc. 
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That  was  our  rather  ambitious  plan  when  we  started  out;  limitations  of  support  and 
available  resources,  particularly  available  authors,  have  determined  how  much,  or  how 
little,  we  could  accomplish.   So  far,  as  indicated  above,  we  have  concentrated  on  Part  I — 
Concepts.   The  Part  I  chapters  already,  or  soon  to  be,  published  are: 

1.  Introduction,  by  F.  E.  Nicodemus,  H.  J.  Kostkowski,  and  A.  T.  Hattenburg 

2.  Distribution  of  Optical  Radiation  with  Respect  to  Position  and 

Direction — Radiance,  by  F.  E.  Nicodemus  and  H.  J.  Kostkowski 

3.  Spectral  Distribution  of  Optical  Radiation,  by  F.  E.  Nicodemus 

and  H.  J.  Kostkowski 

4.  More  on  the  Distribution  of  Optical  Radiation  with  Respect  to 

Position  and  Direction,  by  Fred  E.  Nicodemus 

5.  An  Introduction  to  the  Measurement  Equation,  by  Henry  J.  Kostkowski 

and  Fred  E.  Nicodemus 

6.  Distribution  of  Optical  Radiation  with  Respect  to  Polarization, 

by  John  B.  Shumaker 

7.  The  Relative  Spectral  Responsivity  and  Slit-Scattering  Function 

of  a  Spectroradiometer,  by  Henry  J.  Kostkowski 

8.  Deconvolution,  by  John  B.  Shumaker 

9.  Physically  Defining  Measurement-Beam  Geometry  by  Using  Opaque 

Barriers,  by  Fred  E.  Nicodemus 

10.   Introduction  to  Coherence  in  Radiometry,  by  John  B.  Shumaker. 

Chapters  1,  2,  and  3  were  in  TN  910-1;  4  and  5  in  TN  910-2;  6  in  TN  910-3;  and  7,  8,  and  9 
in  TN  910-4.   In  addition,  the  final  draft  for  chapter  10  (Part  I)  is  now  being  revised 
and  should  be  ready  for  publication  in  TN  910-6  within  one  or  two  months  following  the 
publication  of  this  Technical  Note  910-5. 

Contingent  on  the  availability  of  outside  funding  support,  we  hope  that,  during  the 
next  two  to  three  years,  we  can  continue  the  current  arrangements  with  Catholic  University 
of  America  as  described  above  and  complete  the  following  additional  chapters  of  Part  I — 
Concepts: 

Distribution  of  Optical  Radiation  with  Respect  to  Time,  by  Fred  E.  Nicodemus 

Linearity  Considerations  and  Calibrations,  by  John  B.  Shumaker 

Spectrophotometry,  by  William  H.  Venable,  Jr. 

Blackbody  Radiation  and  Temperature  Scales,  by  Henry  J.  Kostkowski 

Physical  Photometry,  by  A.  T.  Hattenburg 

Thermal  Radiation  Properties  of  Matter,  by  Joseph  C.  Richmond 

Spectroradiometry  of  Spectral  Lines,  by  Henry  J.  Kostkowski. 

As  stated  above,  the  measurement-equation  approach  is  central  to  the  entire  Manual 
and  the  material  presented  in  each  chapter  needs  to  be  related  to  the  approach  introduced 
in  Chapter  5.   Depending  on  how  well  this  is  accomplished  in  each  chapter,  there  may  be 
need  at  the  end  to  have  a  final  summary  chapter  for  this  purpose,  tying  up  loose  ends  and 
putting  the  whole  Part  I  into  perspective.   Also,  examples  of  various  categories  of  envi- 
ronmental and  instrumental  parameters  and  their  significance  could  be  usefully  presented 
and  discussed.   But  it's  too  soon  to  evaluate  this  need  at  present. 


Incidentally,  in  preparing  material  for  the  Manual,  we  have  had  the  pleasure  of 
rediscovering  the  fact  that  the  best  way  to  learn  anything  is  to  try  to  teach  it  to  someone 
else.   The  exercise  of  preparing  tutorial  material  for  such  wide  general  application  has 
required  us  to  analyze  our  own  measurement  activities  in  a  different  way,  broadening  our 
understanding  and  resulting  in  improved  methods  and  more  accurate  results.   Note  that  all 
references  to  measurement  accuracy  or  uncertainty  in  this  preface  are  concerned  not  only 
with  precision  (relating  to  the  repeatability  of  measurement  results)  but  also  with 
accuracy  (relating  to  agreement  with  the  "truth"  which,  while  unknowable  in  the  last 
analysis,  is  approximated  by  analyses  and  estimates  based  on  the  widest  possible  experience, 
including  agreement  with  measurements  of  the  same  quantities  by  others,  particularly  when 
they  have  used  different  instrumentation  and  methods  of  measurement) . 
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SELF-STUDY  MANUAL  on  OPTICAL  RADIATION  MEASUREMENTS 

Part  III.   Applications 

This  is  the  first  chapter  in  Part  III  (Applications)  of 
this  Manual.   Material  developed  in  earlier  chapters  is  used 
to  perform  state-of-the-art  measurements  in  Gainesville,  Florida, 
of  solar  terrestrial  spectral  irradiance  between  290  [nm]  and 
340  [nm] . *   The  measurement  equation  is  used  to  design  the  exper- 
iment and  to  address  the  effects  of  polarization,  non-linearity, 
spectral  scattering,  distortion,  slit-scattering  function,  spec- 
tral irradiance  calibration,  and  wavelength  calibration.   Estimates 
are  made  of  the  uncertainties  associated  with  all  these  factors. 
The  total  uncertainty  is  estimated  to  be  about  10  percent.   An 
appendix  includes  details  for  computing  the  solar  terrestrial  spec- 
tral irradiance  between  280  and  310  [nm] .   Suggestions  are  made  for 
reducing  the  uncertainty  by  about  one  third  and  for  further  research 
in  UV  solar  terrestrial  measurements. 

Key  Words:      atmospheric  attenuation;  atmospheric  ozone;  optical 
radiation  measurements;  radiometry;  solar  radiation;  spectro- 
radiometry;  UV  spectral  measurements. 

Chapter  1.   Measurement  of  Solar  Terrestrial  Spectral  Irradiance  in  the  Ozone 
Cut-Off  Region* 

t 
by  Henry  J.  Kostkowski  ,  Robert  D.  Saunders,  John  F.  Ward, 

Charles  H.  Popenoe,  and  A.E.S.  Green** 

In  this  CHAPTER.   We  use  the  material  developed  in  all  the  previous  chapters  of  this 
Self-Study  Manual  to  illustrate  how  a  very  difficult  radiometric  problem  is  attacked  and 
solved.   The  problem  is  the  measurement  of  solar  terrestrial  spectral  irradiance  between 
about  290  [nm]  and  340  [nm] ,  the  region  where  absorption  by  atmospheric  ozone  produces  a 
sharp  cut-off  in  the  observed  spectrum.   In  addition  to  its  general  tutorial  value,  the 
.chapter  has  applications  in  research  on  the  biological  effects  of  UV  and  in  atmospheric 
ozone  determinations. 

THE  PROBLEM  —  WHAT  and  WHY.   The  optical  radiation  measurement  problem  that  this  chapter 
addresses  is  the  determination  of  the  solar  spectral  irradiance  at  a  point  on  the  earth's 
surface  at  wavelengths  less  than  about   340  [nm] .   A  typical  plot  of  the  spectral  irradi- 
ance (averaged  over  a  wavelength  interval  of   1  [nm])   in  this  region  is  shown  in  figure  1, 
The  sharp  spectral  cut-off  is  caused  by  absorption  in  atmospheric  ozone. 


*Equipment,  experimental  measurements,  and  analyses  associated  with  this  chapter  were 
extensively  supported  by  the  Environmental  Protection  Agency. 

t 
Current  address,  Spectroradiometry  Consulting,  3972  Wilcoxson  Drive,  Fairfax,  Virginia 

22031. 

**Physics  Department,  University  of  Florida,  Gainesville. 

*Unit  symbols  are  all  enclosed  in  square  brackets  to  emphasize  dimensionality  and  encourage 
routine  checks  of  unit-dimension  consistency,  as  discussed  in  Chapter  1  of  Part  I,  on  p.  8 
of  NBS  TN  910-1. 
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Figure  1.   Typical  solar  terrestrial  spectral  irradiance, 
averaged  over  a  wavelength  interval  of  1  [nm] , 
and  the  action  spectrum  for  DNA. 


A  major  application  for  the  measurement  of  solar  terrestrial  spectral  irradiance  in 
the  "ozone  cut-off  region"  occurs  in  the  many  investigations  that  are  being  conducted  on 
the  effect  of  this  radiation  on  living  things.   Optical  radiation  in  this  region  is 
generally  damaging  to  animals  and  plants.   The  action  spectrum  (relative  spectral  effec- 
tiveness) for  destroying  DNA  in  cells  is  shown  in  figure  1.   The  damage  produced  by  an 
additional  increment  of  irradiance  at   310  [nm]   is  a  thousand  times  more  than  at   350 
[nm] .    Therefore  accurate  measurement  of  this  solar  terrestrial  spectral  irradiance  is 
important  in  experiments  investigating  its  biological  effect.   Unfortunately,  accurate 
measurements  of  the  spectral  irradiance  in  this  region  are  very  difficult  because  of  its 
rapid  change  with  wavelength  and  the  very  large  magnitude  of  that  change. 

Another  application  of  solar  terrestrial  measurements  in  this  region  is  in  determining 
the  amount  of  ozone  in  the  atmosphere.   This  quantity  is  usually  obtained  from  measurements 
at  wavelengths  between   305  [nm]   and   340  [nm]   using  the  Dobson  method  [1]  .   Because 
ozone  has  a  much  larger  effect  on  the  shorter  wavelength  radiation,  it  has  been  suggested 
[2,3]  that  measurement  of  the  shape  of  the  solar  terrestrial  spectrum  in  the  "cut-off 
region"  be  used  as  a  means  for  determining  atmospheric  ozone  thicknesses.   A  major  step  in 
pursuing  this  suggestion  is  a  demonstration  that  the  greater  difficulties  inherent  in 
these  measurements  do  not  outweigh  the  advantage  of  increased  sensitivity  for  ozone  deter- 
mination. 

Finally,  the  measurement  of  the  solar  terrestrial  spectral  irradiance  in  the  ozone 
cut-off  region  is  of  interest  to  us  in  this  Manual  because  it  is  a  very  difficult  radio- 
metric problem,  and  it  provides  an  opportunity  to  illustrate  the  application  of  many  of 
the  concepts  and  techniques  introduced  in  Part  I. 

MEASUREMENT  CONFIGURATION.    The  first  step  that  should  be  taken  in  most  optical  radiation 
measurement  problems,  as  discussed  in  Chapter  5,  Part  I — Concepts  (p.  87,  910-2)2  [4],  is 
the  preparation  of  a  diagram  of  the  entire  measurement  configuration.   This  should  show 
the  size  and  location  of  the  unknown  source  and  the  calibration  standard  relative  to  the 
measuring  instrument,  the  location  of  any  sources  of  unwanted  or  background  radiation,  the 
spectral  distribution  and  approximate  magnitude  of  the  radiometric  quantity  of  interest 
for  all  the  sources  (unknown,  standard,  and  background) ,  and  details  of  the  measuring 
instrument  including  size  and  position  of  the  receiving  aperture,  acceptance  solid  angle 
and  spectral  responsivity  function  (Chapters  9  and  7,  910-4)  [5,6]. 

Usually  the  complete  measurement  configuration  will  not  be  known  when  one  begins 
planning  a  measurement  but  will  evolve  as  these  plans  evolve.   The  first  bits  of  informa- 
tion will  often  be  obtained  from  talking  with  colleagues,  from  technical  reports,  or  from 


Figures  in  brackets  indicate  literature  references  listed  at  the  end  of  this  chapter. 

References  to  specific  material  in  various  chapters  of  the  Self-Study  Manual  will  be 
given  in  this  form  —  (page  number  or  chapter,  NBS  Technical  Note  number). 
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the  archival  literature.   In  some  cases,  preliminary  measurements  must  be  made  in  order  to 
get  the  required  information.   The  final  detailed  configuration  is  only  available  after 
the  measurements  have  been  completed.   Even  then,  one  can  often  use  what  was  learned  in 
the  first  measurement  to  change  the  configuration  and  make  a  second  more  accurate  or  more 
efficient  measurement.   We  will  indicate,  in  presenting  our  configuration,  some  of  this 
evolution,  but,  in  general,  the  configuration  described  will  be  the  one  actually  used  in 
the  measurements  we  made  in  Gainesville,  Florida,  during  June  1980. 

Approximate  spectral  irradiance.    The  first  thing  we  sought  relative  to  the  measurement 
configuration  was  an  approximation  to  what  we  were  trying  to  measure  —  the  spectral 
distribution  and  magnitude  of  the  solar  terrestrial  spectral  irradiance.   This  was  obtained 
at  a  solar  symposium  that  we  had  attended  a  few  years  ago  and  is  shown  in  figure  1.   A 
report  by  Bener  [7]  and  a  few  preliminary  measurements  confirmed  the  general  shape  and 
magnitude  shown  in  figure  1.   The  approximation  was  further  refined  by  calculations1  made 
using  data  from  the  literature  on  extraterrestrial  sunlight  in  high  spectral  resolution 
and  typical  transmittances  of  atmospheric  molecules,  particles  (aerosols)  and  ozone. 
Ultimately,  a  set  of  ozone  and  aerosol  data  obtained  in  our  June  measurements  was  factored 
in.   The  distribution  shown  in  figure  2  was  calculated  for  a  spectral  slit  width  of 
1  [nm] .    It  represents  an  average  of  the  computed  spectral  irradiance  weighted  by  a 
triangular  slit-scattering  function  (p.  15,  910-4)  [6]  with  a  (full)  width  at  half-height 
of   1  [nm] .    Using  a  much  smaller  spectral  slit  width   (0.003  [nm]),   the  calculations  of 
the  average  spectral  irradiance  have  the  appearance  shown  in  figure  3.   Figure  4  gives  the 
results  for  spectral  slit-widths  of  0.18  and  0.34  [nm]   with,  in  addition,  two  different 
assumed  conditions  of  atmospheric  attenuation.   In  this  latter  figure  we  have  also  included 
the  spectral  irradiance  of  a  1000  watt  [W]   tungsten  halogen  lamp  standard  at  a  distance 
of   50  [cm]. 

The  important  points  to  note  relative  to  the  solar  spectral  irradiance  distribution 
are: 

1.  The  spectral  irradiance  has  a  great  deal  of  spectral  structure 
(Fraunhofer  lines)  so  that  with  the  spectral  resolution  normally 
available  in  spectroradiometry  the  quantity  determined  will  be 
an  average  weighted  by  the  spectral  responsivity  function  of  the 
measuring  instrument. 

2.  For  a  spectral  slit  width  of   1  [nm]   the  averaging  process 
has  "washed  out"  most  of  the  structure. 

Geometry.  The  next  aspect  of  the  measurement  configuration  we 'examined  was  the  geometry. 
The  set-up  available  to  us  for  making  the  solar  measurements  was  in  a  roof-top  laboratory. 
A  solar  tracker  on  the  roof  directed  the  solar  beam  through  a  hole  in  the  roof  down  to  the 
spectroradiometer,  as  shown  in  figure  5. 


details  of  these  calculations  are  presented  in  Appendix  1A. 
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Figure  2.    Calculated  solar  terrestrial  spectral  irradiance  (triangular 
slit-scattering  function  with  a  1  [nm]   spectral  slit  width) 
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Figure  3.    Calculated  solar  terrestrial  spectral  irradiance 
(0.003  [nm]   spectral  slit  width). 
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Figure  4.    Calculated  solar  terrestrial  spectral  irradiance   (0.18  and  0.34  [nm] 
spectral  slit  widths)  and  a   1000  [W]   tungsten-halogen-standard-lamp 
spectrum  at   50  [cm].    (Vertical  separation  of  the  solar  curves  is 
due  to  different  assumed  atmospheric  attenuation  conditions.) 
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The  tracker  consists  of  two  plane  mirrors  (see  figure  15  as  well  as  figure  5) .   The 
two  mirrors  can  be  rotated  as  a  unit  about  a  vertical  axis  through  the  center  of  the 
mirror  closest  to  the  monochromator,  and  the  mirror  farthest  from  the  monochromator  can  be 
rotated  about  a  horizontal  axis  through  the  centers  of  both  mirrors.   Both  rotations  are 
computer  controlled  [3].   The  various  parts  of  the  spectroradiometer  shown  in  the  block 
diagram  of  figure  5  will  be  discussed  later. 

Since  we  were  interested  in  determining  the  spectral  irradiance  of  the  direct  sun 
rather  than  the  global  irradiance  (sun  plus  sky),  we  had  to  limit  our  field  of  view  so  that 
we  included  negligible  sky  irradiance.   In  order  to  minimize  the  effect  of  any  tracking 
errors,  we  chose  a  3.25-inch  [in]   aperture  resulting  in  a  sufficiently  large  field  of 
view  (see  figure  6)  relative  to  tracking  errors  but  one  for  which  we  thought  the  sky 
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Figure  6.   Measurement  configuration  (without  a  diffuser) 


irradiance  would  be  negligible,   1.5  [deg].   Measurements  made  by  blocking  out  the  sun 
with  opaque  discs  of  various  diameters  confirmed  that  the  sky  irradiance  was  negligible1 
with  the  1.5  [deg]   field  —  about  0.5  percent  that  of  the  sun. 

Measurement  with  and  without  a  diffuser.    Normally,  spectral  irradiance  measurements  are 
made  with  a  diffuser  (e.g.  an  averaging  sphere)  as  the  first  optical  element  in  the 
spectroradiometer  (p.  71,  910-2)  [4],   The  purpose  of  this  is  to  insure  that  the  spectral 
flux  responsivity  is  invariant  with  direction  of  the  incident  flux.   In  this  case,  the 
measurement  equation  (p.  65,  910-2)  [4] 


s(xo)  = 


R  •!■  •cos6«du)«dA«dX  (1.1) 

i  ■    <P   A 
AAJ  AJu 


may  be  reduced  to 


s(xo)  = 


R  -E  «dA-dX,  (1.2) 

AAJA 


which  is  the  equation  required  for  a  measurement  of  spectral  irradiance  E. .   We  used  a 

A 

1/32  [in]   thick  piece  of  teflon  as  a  transmitting  diffuser,  positioning  it  so  that  it 
approximately  filled  the  field  of  view  of  the  monochromator.   The  geometry  of  our  measure- 
ment configuration  (not  to  scale)  using  this  diffuser  is  shown  in  figure  7. 

There  is  one  enormous  disadvantage  of  the  diffuser.   It  reduces  the  flux  reaching  the 
detector  by  a  large  factor.   In  our  case  this  was  about   1000.   Therefore,  below  305 
[nm] ,   where  we  were  beginning  to  be  energy  limited,  we  rarely  used  a  diffuser.   The 
entrance  slit  and  the  sun  itself  defined  the  beam  geometry  as  shown  in  figure  6.   In  this 
case,  we  had  to  rely  on  the  flux  responsivity  being  sufficiently  uniform  over  the  0.5 
[deg]   field  of  view  without  a  diffuser.   It  was  essential  to  confirm  that  this  was  the 
case  by  measuring  the  flux  responsivity  as  a  function  of  angle.   This  was  done  by  using  a 
small  tungsten  halogen  lamp  providing  a  beam  entering  the  monochromator  with  an  angular 
extent  of   0.2  [deg]   in  one  direction  and  of   0.6  [deg]   perpendicular  to  this  direction. 
R  (0)   was  obtained  by  determining  the  output  signal  as  a  function  of  angle  of  tilt  of  the 
spectroradiometer.   The  result,  shown  in  figure  8,  was  about  the  same  for  either  orienta- 
tion of  the  lamp.   The  responsivity  was  uniform  around  the  peak  to  within  1  percent 
between  about  -0.25  [deg]   and  +0.25  [deg].    This  meant  that  the  standard  should  be 
confined  to  a  0.5  [deg]   field  of  view,  the  angle  subtended  at  the  earth's  surface  by  the 
sun.   This  was  accomplished  by  positioning  the  approximately  14  [mm]  x  36  [mm]   standard 


Negligible  is  in  reference  to  the  overall  or  total  uncertainty  that  we  were  striving  to 
realize,  which  in  these  very  difficult  solar  spectral  measurements  was  about   5  or  10 
percent. 
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APERTURE 


TEFLON  DIFFUSER 


Figure  7.    Measurement  configuration  with  a  diffuser, 


lamp  on  the  solar  side  of  the  tracker  about   400  [cm]   from  the  entrance  slit  of  the 

spectroradiometer.   It  was  also  necessary  that  we  be  able  to  align  and  maintain  the 

spectroradiometer  orientation  relative  to  the  beam  from  the  tracker  to  within   0.25  [deg] . 

For  this  purpose  a  mount  was  designed  for  the  spectroradiometer  which  provided  micrometer- 

_3 
type  adjustments.   The  precision  of  the  angular  positioning  was   10  ~  [rad]   (0.05  [deg]) 

and  the  translational  positioning   0.1  [mm]. 

In  addition  to  the  above  described  check  on  the  uniformity  of  the  flux  responsivity, 
we  decided  to  confirm  by  direct  measurement  that  not  using  a  diffuser  did  not  greatly 
degrade  the  accuracy.   The  confirmation  could  be  made  in  the  wavelength  region  between 
305  [nm]  and  340  [nm]   where  the  spectroradiometer  output  signal  was  high  even  when  the 
diffuser  was  used.   We  planned  to  make  several  sets  of  measurements  in  this  region  with 
and  without  the  diffuser.   If  significant  differences  were  found,  we  could  then  consider 
using  these  differences  to  correct  the  shorter  wavelength  data  where  the  signal  was  too 
small  to  use  the  diffuser. 
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Figure  8.   Angular  responsivity  function  of  NBS  Solar 
Spectroradiometer  (without  dif fuser) . 
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Calibration  standard.    The  spectral  irradiance  source  standard  usually  used  to  calibrate 
spectroradiometers  in  the  wavelength  region  in  which  we  were  working   (287-340  [nm])   is  a 
1000  [W]   clear  tungsten  halogen  lamp.   Figure  4  includes  a  wavelength  plot  of  the  spectral 
irradiance  at   50  [cm]   from  such  a  standard.   The  spectral  irradiance  of  this  standard  is 
quite  different  from  that  of  the  sun  for  most  of  the  wavelength  range  of  interest.   Also 
the  beam  geometry  (coiled-coil  filament  is  about   5  [mm]  x  25  [mm])   is  different  from 
that  of  the  sun.   This  (standard)  flux  also  has  a  polarization  of  about  3  percent  while  the 
sun  is  unpolarized.   Differences  between  the  optical  radiation  being  measured  and  that  of 
the  standard  are  potential  sources  of  error  [8] .   In  characterizing  the  measuring  instru- 
ment, which  will  be  discussed  later,  the  effects  of  such  differences  will  have  to  be 
determined. 

MEASUREMENT  EQUATION.    The  general  measurement  equation  for  incoherent  radiation  including 
polarization  effects  (p.  23,  910-3)  [9]  is 


S(Aq)   =  (Roo'L0  +  R01*L1  +  R02*L2  +  R03*VCOse*da)'dA'dX'      (1,3) 

AX-*  AA-J  Aw 


A  similar  equation  is  required  for  the  standard  which  in  our  case  would  be  located  on  the 
solar  side  of  the  aperture  shown  in  figure  5.   This  equation  is 


S-(Ao)      = 


AA 


AAJ  Aoo 


s    (R00'L0  +  R01-Ll  +  R02*L2S  +  R03-L3S)-COSe'd^dA-dA 


(1.4) 


where  the  R_. .   are  the  elements  in  the  first  row  of  a  Mueller  responsivity  matrix  repre- 
Oi  J 

senting  a  combination  of  the  tracker  and  spectroradiometer  responsivities;  that  is  (Chap- 
ter 6,  910-3)  [9], 
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(1.5) 


where  the  L.   are  the  components  of  the  Stokes  vector  associated  with  the  sun  or  the  stand- 

g 
ard,   L  .    Since  the  sun  is  unpolarized  L..  =  L„  =  L„  =  0   and 


s(xo)  = 


AAJ AA;  Ao) 


Rn»L  «cos6  'dwdA-dA. 


(1.6) 


Clearly,  if  the  usual  3  percent  polarized  standard  lamp  were  used,  we  would  have  to  concern 
ourselves  with  a  more  complex  measurement  equation  for  the  standard,  one  involving  other 
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R„  .   and  L..    Therefore  we  chose  to  calibrate  a  frosted  (unpolarized)  tungsten  halogen 
lamp  by  comparing  it  to  a  typical  standard,  using  a  (different)  polarization- insensitive 

spectroradiometer,  and  use  this  unpolarized  secondary  standard  in  our  solar  measurements. 

s     s     s 
The  measurement  equation  for  such  an  unpolarized   (L1   =  L„   =  L„   =0)   standard  would 

then  also  be 


ss(xo)  = 

'AAJAA 


R_n.Ln«cos6«da)»dA'dX.  (1.7) 

,  s  00  0 

Au) 


Thus  it  appears  that  since  both  the  sun  and  the  standard  lamp  are  unpolarized,  there  are 
no  effects  due  to  polarization.   However,  this  is  not  the  case.   The  first  column  of  the 
responsivity  matrix  for  the  tracker  in  eq.  (1.5)  is  a  function  of  the  tracker  orientation 
[10]  resulting  in  R„n  also  changing  with  tracker  orientation.   This  orientation  is 
determined  by  the  time  of  day  that  the  solar  measurements  are  made.   Therefore  the  (com- 
bined) measurement  equation  should  be  written 

S(X   )  !         Ron(T)«Ln-cos6«dwdA-dX 

°  _      •lAXJAA-lAa)  ,,    QN 

SS(A    )  iff  RAA(T).L^-cose-dco.dA.dA 

°  L  L  L    s     00  0 

1  A\J  &AJ  Aa) 

and  the  calibration  must  be  obtained  as  a  function  of  tracker  orientation,  expressed  here 
in  terms  of  the  time   (T)   at  which  the  instrument  correctly  points  to  the  sun,  so  that 
the  R„n(T)   corresponding  to  the  time  of  the  solar  measurement  is  used.   Figure  9  shows 
how  R„n  varies  as  a  function  of  tracker  orientation  from  290  to  300  [nm] .    From  300 
to  340  [nm]   the  functional  relationship  is  different  than  that  at  the  shorter  wavelengths 
and  also  varies  with  wavelength. 

When  the  diffuser  is  used,  calibrations  need  not  be  taken  as  a  function  of  tracker 
orientation.   In  this  case  the  spectroradiometer  is  unpolarized  (all  R. .   in  the  spectro- 
radiometer 's  responsivity  matrix,   except  ^nn,      are  equal  to  zero),  and  since  the 
R„n(tracker)   varies  only  slightly   (0.1  percent)   with  tracker  orientation  [10]  so  does 
R   (spectroradiometer  and  tracker).    In  other  words,  although  the  tracker  polarizes  the 
incoming  radiation,  it  doesn't  change  the  flux  level  significantly.   Therefore  when  the 
spectroradiometer  is  polarization  insensitive  (when  using  the  diffuser) ,  the  responsivity 
is  virtually  independent  of  tracker  orientation. 

Since  we  wish  to  measure  spectral  irradiance,  the  measurement  eq.  (1.8)  must  be 
modified  to  include  E   rather  than  L  .   As  discussed  earlier,  this  requires  that 

A  A 

g 

Rnn(T)   be  independent  of  the  source  solid  angles   Aco  and  Aw  .    Then 
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Figure  9.    Relative  system  responsivity  (without  diffuser)  as  a 
function  of  tracker  mirror  orientation   (T) .   (Time 
of  day  T   is  a  convenient  variable  for  designating 
the  orientation  of  a  mirror  that  is  tracking  the  sun.) 
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s(xo) 


ss(xo) 


R00(T)-(    L0-cos6-dco)-dA-dX 
AXJAA        ^Ao) 


AX 


R00(T)«M    LQ»cos6»daj)  «dA«dX 

AA        V-'Acj 


(1.9) 


and 


s(xo) 


ss(xo) 


AX 


AA 


Roo(T)'VdA*dx 


Roo(T)-Ex*dA'dx 


(1.10) 


AXJAA 


In  our  application,   E  and  E   are  constant  over  the  receiving  aperture  area  AA  to 
well  within  one  percent  so  that 


s(xo) 


ss(xo) 


R00(T).dA].dX 

AX    \JAA 


R00(T).dA).dX 

AX    \JAA 


(1.11) 


s(xo) 


dX 


Lwt)- 

SS(Xo)      f  E^-Rp(T).dX 

Jax  x    e 


(1.12) 


-1  2 
where  Rp(T)  [S*W   «m  ]   is  the  irradiance  responsivity  expressed  as  a  function  of  tracker 

orientation,   T.    Under  certain  conditions,  which  will  be  discussed  later, 


S(Xq)      R^(Xo,T).[   EA-z(Xo,X,T)-dX 

•'ax 


(1.13) 


S(V       ^W' 


AX 


E  -z(X  ,X,T  )-dX 
X    o    s 


where  R_(X  ,T)  [S«W   'in  ]   is  the  irradiance  responsivity  factor1  at  wavelength  X   and 


*In  NBS  Tech.  Note  910-4  [6]  (pp.  15  and  17),  the  responsivity  factor  r   was  defined  in 
terms  of  the  relative  responsivity  r  =  z«r  ,   and  R^   =  K«r  =  K«z«r  .    In  this  chapter 
it  is  more  appropriate  to  absorb  K  into  the  responsivity  factor  resulting  in  R^  =  z«R^, 
where  R^   is  called  the  irradiance  responsivity  factor.   Note  that  R^,   K,   and  R^ 
all  have  the  same  unit-dimensions   [S«W  «m  ],   while  r,   r  ,   and  z  are  dimensionless. 
In  general,  the  slit-scattering  function  is  also  a  function  of  the  polarization  of  the 
incident  flux  and  therefore  of  the  tracker  orientation  T,   as  indicated  in  eq.  (1.13). 
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tracker  orientation  T,  and  where  z  is  the  slit-scattering  function  (p.  15,  910-4)  [6]. 
When  the  diffuser  is  employed,  R^  or  It  and  z  will  be  independent  of  tracker  orien- 
tation T. 

OTHER  FACTORS  in  EXPERIMENTAL  DESIGN 

Scattering  and  distortion.    For  many  spectroradiometric  measurements,  the  error  produced 

by  assuming  that  E..   is  constant  (so  that  it  may  be  moved  outside  of  the  integral  sign) 

A 

in  eq.  (1.12)  or  eq.  (1.13)  is  not  great.   This  is  not  so  in  our  case. 

The  biggest  problem  is  the  unwanted  response  to  the  large  and  extensive  (relative  to 
wavelength)  spectral  flux  far  from  the  wavelength  setting  of  the  spectroradiometer  that  is 
detected  by  the  instrument  (p.  24,  910-4)  [6],   This  is  caused  by  the  very  long  wavelength 

wing  of  the  responsivity  function  (see  figure  10) .   If  this  wing  of  the  responsivity 

—8 
function  is  even  10    relative  to  its  peak,  a  typical  value  for  a  high-quality  double 

monochromator,  calculations  show  that  the  amount  of  flux  "leaked"  or  "scattered"  at  290 
[nm]   for  a  spectral  distribution  similar  to  that  in  figure  1  is  about  equal  to  the  actual 
spectral  irradiance  at  290  [nm] .    The  amount  of  leakage  changes  very  little  with  mono- 
chromator wavelength  setting.   Thus  we  expected  significant  scattering  below  about  293 
[nm]   and  would  probably  have  to  determine  the  wings  of  the  responsivity  function  exper- 
imentally and  deconvolute  (p.  35,  910-4)  [11]  the  numerator  of  eq.  (1.12)  if  results 
accurate  to  5-10  percent  are  desired  at  the  short  wavelengths. 

In  addition  to  the  wings  of  the  responsivity  function  producing  an  error  or  uncer- 
tainty due  to  scattering,  the  central  portion  can  cause  an  error  if  the  spectral  irradiance 
being  measured  does  not  vary  linearly  with  wavelength  (p.  26,  910-4)  [6].   This  effect  is 
referred  to  as  distortion.   Calculations  show  that  the  curvature1  of  the  spectral  distribu- 
tion in  figure  2  is  sufficient  to  produce  an  error  of  about  20  percent  when  the  spectral 
slitwidth  is  1  [nm] .   However,  if  one  were  only  interested  in  the  average  spectral  irra- 
diance over  a  1  [nm]   interval,  the  error  would  be  reduced  to  about  5  percent.   Because 
•of  the  structure  actually  present  in  the  spectral  distribution  (figures  3  and  4) ,  use  of 
a  narrower  spectral  slit  width  such  as  0.18  [nm]   does  not  reduce  the  distortion.   There- 
fore, deconvolution  calculations  (Chapter  8,  910-4)  [11]  of  some  type  will  be  required  and 
for  this  purpose  the  central  portion  of  the  responsivity  function  will  have  to  be  deter- 
mined experimentally. 

Wavelength  accuracy.    A  disadvantage  of  using  a  1  [nm]   spectral  slit  width  as  in  figure 
2  is  that  the  wavelength  uncertainty  of  the  instrument  would  have  to  be  limited  to  about 
0.01  [nm] .    The  slope  of  the  spectral  distribution  in  figure  2  at  290  [nm]   is  about  100 
percent  per  nanometer  [nm] .    The  wavelength  uncertainty  of  the  best  small-  or  medium- 
sized  monochromators  is  usually  10  to  20  times  larger  than  0.01  [nm] .    It  appears  that 
obtaining  adequate  wavelength  accuracy  will  be  a  problem  in  this  measurement.   One 
solution  would  be  to  use  a  small  spectral  slit  width  and  use  the  features  of  the  spectral 


*To  see  the  extensive  curvature  one  must  replot  the  figure  with  a  linear  ordinate. 
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Figure  10.   Response  in  the  cut-off  region  showing  the  effect  of 
the  long  wavelength  wing  of  the  responsivity  function. 
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structure  seen  in  figure  4  as  a  (wavelength)  self-calibration  to  help  reduce  the  wave- 
length uncertainty. 

Computer  control.    Finally,  careful  analysis  of  the  measurement  configuration  has  shown 
that  without  a  diffuser  —  and  this  is  the  configuration  required  for  the  low  spectral 
fluxes  below   305  [nm]  —  the  responsivity  function  will  depend  on  time  (through  tracker 
orientation)  as  indicated  in  eq.  (1.12).   In  addition,  temporal  variations  in  atmospheric 
conditions  (e.g.,  clouds)  make  it  desirable  to  perform  the  measurements  in  as  short  a 
time  as  possible.   On  the  other  hand,  the  signal  at  the  shortest  wavelengths  will  be  very 
small  and  probably  require  integrating  for   5  or  10  [min] .    The  flexibility  required  to 
handle  these  various  timing  problems  efficiently  suggests  utilizing  a  computer  for  con- 
trolling the  time  sequence  in  which  data  are  taken  and  for  storing  and  conveniently 
reducing  these  data. 

DESCRIPTION  of  SPECTRORADIOMETER.    After  careful  consideration  of  the  various  factors 
discussed  in  the  previous  sections,  the  spectroradiometer  system  shown  in  figure  5  was 
designed  and  developed.   There  are  four  major  subsystems  of  the  spectroradiometer  con- 
sisting of  the  monochromator ,  the  detection  system,  the  amplification  and  counting  system, 
and  the  computer  control  and  data  handling  system. 

The  major  monochromator  specifications  that  were  either  required  or  highly  desirable 
were: 

(1)  Wavelength  reproducibility  of   0.01  [nm] . 

(2)  Spectral  slit  width  of  about   0.2  [nm]   and  possibly   1  [nm] . 

—8 

(3)  A  slit-scattering  function  that  did  not  exceed  10 

of  its  peak  value  at   50  [nm]   from  its  peak  wavelength. 

(4)  Angular  aperture  at  least  as  large  as  that  of  an  f/5  system. 

(5)  Spectral  range  at  least  from  280  [nm]   to  370  [nm] . 

(6)  Weight  less  than   1  [kg]   and  dimensions  less  than 
25  [cm]  x  25  [cm]  x  50  [cm]. 

We  were  fortunate  to  find  that  Dr.  William  Fastie  of  Johns  Hopkins  University  had 
recently  designed  a  monochromator  intended  to  meet  or  exceed  these  specifications.   A  few 
of  these  instruments  had  already  been  produced  commercially.1   We  borrowed  one  of  these 
monochromators  and  checked  its  performance  relative  to  the  above  specifications.   All 
specifications  were  satisfied  except  the  wavelength  reproducibility  which  was  about  five 
times  poorer  than  desired.   We  made  some  design  modifications  that  we  thought  had  a  good 
chance  of  improving  the  wavelength  performance  and  let  a  contract  for  construction.   The 
major  features  of  the  final  design  were: 

(1)  All  Invar  36  construction  and  all  optical  materials  of  Cervit 
to  assure  thermal  stability. 

(2)  Screw  driven  sine  wavelength  drive  with  stepping  motor  and 
encoder  readout  to   0.001  [nm]   and  setability  to   0.005  [nm] . 


■^Research  Support  Instruments,  Inc.,  10610  Beaver  Dam  Road,  Cockeysville,  MD  21030. 
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(3)  Double  Ebert-Fastie  (dispersion  cancelling)  monochromator  with 
focal  length  of  approximately  125  [mm]   and  optical  speed  of   f/5, 

(4)  Holographic  gratings  of   3600  [grooves  per  nm] . 

(5)  All  optical  elements  super  polished  and  mirrors  coated  in 
such  a  way  as  to  optimize  high  specular  reflectance  and 
low  diffuse  reflectance  (scattering). 

(6)  Fixed  (non- adjustable)  slits  with  spectral  slit  widths 
of  0.2  [nm]   and   1  [nm] . 

(7)  Weight  less  than  4.5  [kg]   and  dimensions  less  than 
15  [cm]  x   15  [cm]  x  30  [cm]. 

The  cost  was  approximately  $27,000  in  1980. 


RETAINING   RING 
&   ENTRANCE   APERTURE 
WHEN   DIFFUSER   IS   USED 
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Figure  11.    Schematic  diagram  of  dif fuser/f ilter  system. 


20 


The  above  monochromator  did  not  include  a  means  for  mounting  diffusers  or  filters  in 
front  of  the  entrance  slit  or  a  cooled  detector  following  the  exit  slit.  These  additions 
were  designed  and  fabricated  at  NBS. 

The  dif fuser/f ilter  mounting  system  allowed  for  rapid,  reproducible  placement  or 
interchange  of  filters  and  provided  a  means  of  easily  locating  the  optic  axis  of  the 
instrument  for  alignment  purposes.   Figure  11  is  a  schematic  diagram  of  the  system.   The 
tube,  shown  in  the  figure,  is  mounted  by  screws  onto  the  front  (entrance-slit  portion)  of 
the  monochromator  in  a  position  so  that  it  is  coaxial  with  the  monochromator 's  optic 
axis.1   The  cap,  when  screwed  tightly  on  the  tube,  conveniently  insures  good  reproduci- 
bility in  positioning  filters.   One  normally  would  have  as  many  caps  as  filters  or  combi- 
nations of  filters. 

A  schematic  diagram  of  the  detection  system  is  shown  in  figure  12.   The  quartz  lens 
and  photomultiplier  are  positioned  so  that  the  exit  slit  is  defocused  on  the  cathode  of 
the  photomultiplier.   The  larger  area  of  irradiation  minimizes  the  effects  (noise  and 
instability)  of  cathode  non-uniformity  and  of  small  vibrations  or  motions. 


QUARTZ  LENS 


THERMOELECTRIC  COOLER 


MAGNETS 
QUARTZ  WINDOWS 


Figure  12.   Detector  system. 


^he  optic  axis  is  the  central  ray  of  the  full  optical  beam  that  can  traverse  the  mono- 
chromator.  Normally  for  a  well-aligned  monochromator,  the  optic  axis  passes  through  the 
centers  of  all  the  slits  and  the  centers  of  all  the  optical  elements.   A  small  HeNe  laser 
is  useful  for  establishing  the  optic  axis  and  the  alignment  of  this  axis  relative  to  the 
tracker  or  external  apertures. 
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The  purpose  of  the  ring  magnets  is  to  reduce  the  effective  size  of  the  cathode.   The 
magnetic  field  produced  by  the  rings  permits  only  those  thermal  electrons  that  are  emitted 
from  about  a   0.5  [in]   diameter  area  of  the  cathode  to  strike  the  first  dynode.   This 
factor  plus  cooling  to  about   -20°C  with  the  thermoelectric  cooler  minimizes  dark  current 
(or  dark  counts)  which  is  particularly  important  at  the  low  flux  levels  at  wavelengths 
below  about   292  [nm] . 

The  photomultiplier  has  a  bialkali  cathode  and  was  selected  for  high  quantum  effi- 
ciency and  low  dark  count.   In  order  to  reduce  the  amount  of  scattered  signal  detected,  a 
Cs-Te  (solar  blind)  cathode  was  considered.   However,  a  bialkali  cathode  together  with  a 
"solar  blind"  filter  (transmittance  shown  in  figure  13)  accomplishes  the  same  thing  but 
with  a  higher  quantum  efficiency  at   290  [nm] . 

The  output  signal  of  a  photomultiplier,  resulting  from  an  incident  photon  that  has 
ejected  a  photo-electron,  consists  of  a  burst  or  pulse  of  electrons  lasting  only  a  few 
nanoseconds.   When  the  photon  rate  (incident  flux  level)  is  high  these  occur  too  close 
together  (may  even  be  superimposed)  for  them  to  be  detected  individually.   The  average 
electron  flow  must  be  measured.   This  is  the  so-called  current  measuring  mode  which  is 
used  in  most  radiometric  applications  using  photomultipliers.   However,  when  the  flux 
level  is  low,  as  in  our  case  at  the  shorter  wavelengths,  these  electron  bursts  or  pulses 
are  sufficiently  separated  so  that  they  can  be  detected  individually.   This  is  the  photon 
counting  mode.   Because  the  flux  level  at  wavelengths  below  295  [nm]   is  so  low,  we 
decided  to  use  the  photomultiplier  in  a  photon  counting  rather  than  a  current  measuring 
mode. 

The  photon  counting  mode  has  a  number  of  advantages  over  the  current  measuring  mode. 
Noise  and  dark  counts  (signal  counts  when  a  shutter  is  in  place  so  that  the  flux  being 
measured  is  blocked  or  shuttered)  can  usually  be  reduced,  relative  to  the  current  mode, 
because  most  pulses  that  are  produced  by  sources  other  than  incident  photons  (e.g. 
thermionic  electrons  and  stray  electrons  from  dynodes)  have  a  smaller  amplitude  and  can 
be  rejected  by  proper  setting  of  the  threshold  voltage  in  the  discriminator  (see  figure 
5).   Noise  from  electronic  components  and  connections,  leakage  currents,  power  lines,  and 
external  fields  can  be  made  insignificant  or  minimized  through  digital  processing  and 
counting  procedures.   The  average  signal  in  the  current  measuring  mode  is  directly  affected 
by  all  the  above  types  of  noise  but  the  number  of  pulses  per  second  in  the  photon  counting 
mode  is  not.   Another  advantage  in  the  photon  counting  mode  is  that  the  counting  rate  is 
not  as  highly  sensitive  to  the  dynode  voltages  as  is  the  average  signal  in  the  current 
measuring  mode.   Thus,  much  less  stability  and  control  are  required  in  the  high  voltage 
power  supply  used  for  supplying  the  dynode  voltages. 

One  disadvantage  of  the  photon  counting  mode  is  that,  for  very  high  photon  fluxes, 
the  observed  counting  rate   S'  [cs   ]   is  not  exactly  proportional  to  the  incident  flux. 
This  is  because  there  is  a  significant  probability  of  two  photons  arriving  within  the 
resolving  time  (dead  time)  of  the  instrument  (pulse  overlapping) .   Without  this 
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Figure  13.    Spectral  transmittance  of  solar  blind  filter, 
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proportionality  (linearity),  we  can't  do  many  of  the  things  with  the  measurement  equation 
that  we  would  like  to  do.1 

Fortunately  there  is  an  easy  way  of  handling  the  non-linearity  in  the  current  problem. 
This  is  done  by  assuming  that  the  flux  responsivity  is  constant  (linearity  exists)  and 
then  correcting  the  observed  output  signal  to  what  it  would  have  been  if  this  constancy 
had  actually  existed.   Operationally  one  might  think  of  the  "detector"  as  consisting  of 
the  photomultiplier  and  associated  electronics  together  with  a  computer  or  hand  calculator 
which  provides  a  correction  from  the  indicated  deficient  rate  S'   to  a  quantity  S   that 
is  exactly  proportional  to  the  actual  incident  photon  flux.   An  approximate  relationship 
between  S'   and  S   is  [12] 

S/S'   =  C(S')  =  1  +  t-S'  (1.14) 

where  C(S')   is  the  "correction  factor"  for  transforming  the  observed  counting  rate 

S'  [c«s   ]   to  the  "correct"  counting  rate   S,   and  t  [s*c   ]   is  the  resolving  time  of 

the  system.   Differences  between  S'   and  S  can  be  quite  significant.   For  example,  with 

—  7      —1  ft  —  1 

t  =  10    [s«c   ],   a  typical  value,  and  S'  =  10   [cs   ],   the  "correct"  counting  rate 

differs  from  the  observed  counting  rate  by  10  percent. 

High  quality  photon  counting  systems  are  now  available  commercially  as  "shelf  items." 
We  selected  one  of  these  for  this  application.   The  major  points  to  consider  in  this 
selection  are  over-all  stability  and  ease  of  operation  and  required  resolving  time  and 
integration  time. 

The  computer  control  and  data  handling  system  is  the  remaining  subsystem  (figure  5) 
to  be  discussed.   This  subsystem  consists  of  a  48K  RAM  microcomputer,  dual  floppy 
8  [in]   disk  drive,  terminal/video  screen  and  printer.   We  purchased  standard,  readily 
available  commercial  items.   The  electronic  circuitry  (interfacing)  required  to  couple 
the  computer  system  to  the  rest  of  the  spectroradiometer  was  designed  and  constructed  at 
NBS.   This  is  not  a  difficult  task  for  a  digital  electronic  engineer.   It  is  also  possible 
to  obtain  the  interface  electronics  commercially. 

The  computer  can  be  programmed  to  start  measurements  at  a  particular  wavelength, 
integrate  the  counts  for  any  time  interval  from  a  fraction  of  a  second  to  hours  and 
repeat  this  at  any  other  wavelength  or  group  of  wavelengths.   In  addition  to  the 


*When  the  measurement  equation  was  first  derived  in  Chapter  5  [4],  we  indicated  that  the 
assumed  constancy  of  the  flux  responsivity  with  respect  to  the  magnitude  of  the  flux  was  a 
convenience.   However,  we  have  come  to  realize  that  it  is  a  necessary  assumption  for  much 
that  we  want  to  do  when  applying  the  measurement  equation.   Examples  of  this  are  in  going 
from  eq.  (1.12)  to  eq.  (1.13)  and  in  the  definition  R  =  K>z-r  .   (See  footnote  on  page 
16.)   An  in-depth  treatment  of  linearity  is  planned  for  a  forthcoming  chapter  of  this 
Manual . 

24 


wavelength,  total  signal  counts,  and  counts  per  second,  the  ambient  temperature,  time, 
and  date  are  also  recorded.   As  data  are  obtained  they  may  be  viewed  on  the  video  screen, 
and/or  printed  and  stored  on  a  floppy  disk.   Printing  the  data  is  usually  delayed  until 
all  observations  for  the  day  are  complete  because  concurrent  printing  increases  the  time 
required  to  perform  the  experiment. 

CHARACTERIZATION  of  SPECTRO RADIOMETER.    In  order  to  minimize  the  errors  or  uncertainties 
of  a  radiometric  measurement-,  the  measuring  instrument's  responsivity  should  be  determined 
as  a  function  of  each  radiation  parameter  for  which  there  are  significant  differences 
between  the  unknown  radiant  flux  being  measured  and  that  of  the  standard.   The  pertinent 
radiation  parameters  were  discussed  in  the  sections  on  the  measurement  configuration  and 
the  measurement  equation.  The  two  radiation  fluxes  in  this  problem  were  the  sun  and  a 
frosted  1000  [W]   tungsten  halogen  lamp  at  400  [cm];   the  parameters  for  which  the  two 
fields  differed  significantly  were  magnitude  (absolute  value  of  the  spectral  irradiance) 
and  spectral  distribution. 

Relative  to  magnitude.    As  pointed  out  earlier,  the  responsivity  of  photomultipliers 
when  used  as  photon  counters  decreases  at  sufficiently  high  levels  of  incident  photon 
flux  in  accordance  with  eq.  (1.14).   With  the  very  large  magnitude  and  variation  of 
magnitude  of  the  solar  terrestrial  irradiance  in  the  spectral  region  of  interest,  these 
decreases  (non-linearity)  could  be  significant.   This  problem  has  been  addressed  by 
determining  C(S')   of  eq.  (1.14)  and  using  the  "corrected"  output  signal   S  =  C(S')«S' 
instead  of  the  observed  signal  S'   in  the  measurement  equation.   In  this  way  linearity 
with  all  its  desirable  features  is  effectively  maintained. 

The  correction  C(S')   was  obtained  by  first  determining  the  spectral  transmittance 
t  of  a  filter  at  a  sufficiently  low  incident  flux  cj> .   (and  corresponding  low  observed 
counting  rate  S')   where  C(S')  =  1  and  then  repeating  the  measurement  of   t  with 
higher  incident  flux  levels  up  to  the  highest  flux  level  or  counting  rate  obtained  in  the 
solar  measurements.1   In  order  to  achieve  a  more  accurate  approximation,  the  power  series 
in  S'   [eq.  (1.14)]  was  carried  to  one  more  term;  i.e., 

C(S')  =      S/S'   =   1  +  t-S'  +  u-S'2.  (1.15a) 


^-This  technique  can  be  used  in  the  current  problem  because  it  is  well  established  that 
for  photomultipliers,  when  the  counting  rate  is  sufficiently  low,   S'  =  R-$  and  therefore 

$    St  n 

t  =  —  =  -r-f-.   Otherwise,  if  the  responsivity  of  a  detector  were  such  that  S'  =  R«$   at 

i    i 

low  levels,   S'/S'.   would  still  be  constant  as  a  function  of   S!   at  sufficiently  low 
-  t  i  i 

q  »  „ 

i  but  equal  to   x    instead  of   x.   These  filter  measurements  alone  will  not  determine 
the  value  of  the  exponent  n. 
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The  transmittance  t  of  the  filter  referred  to  above  would  then  be 

t   =  C'(SJ)  •  S^/[c(Sp  •  Sj]  (1.16a) 

For  S'   sufficiently  small  C(S!)  =  C(S')  =  1  and  t  =  S'/S^.   For  larger  values  of  S|, 
substituting  eq.  (1.15a)  into  eq.  (1.16a) 

t  =  (s't  +   fS'2  +  u-S^3)/(s^  +  t-S^2  +  u-S^3)  (1.16b) 

and  rearranging  terms 

tS^  -  Sj.  =   (S^2  -  T-S^2).t  +  (S^3  -  x-S^3)-u.  (1.16c) 

Measurements  of  t,   S!   and  S'   were  performed  with  the  253.7  [nm]   spectral  line 
from  a  temperature-stabilized  low-pressure  mercury-discharge  lamp  source  from  which  the 
irradiance  (flux)  level  was  varied  by  the  insertion  of  additional  filters,  taking  care 
not  to  disturb  the  geometry  of  the  beam  incident  on  either  the  filter  being  measured  or 
on  the  monochromator.  Using  these  data,   t  and  u  were  obtained  from  a  least  squares 
fit  of  eq.  (1.16c).   For  values  of  S'   as  high  as  2.5*10   [c«s~  ], 

S/S'   =  C(S')  =  1  +  6.927«10~8«S'  +  7.174-10_15-S'2        (1.15b) 

where  S  is  the  "corrected"  output  signal  and  S'   is  the  actually  observed  output 
signal  in   [cs   ]. 

There  was  another  problem  relative  to  the  responsivity  that  we  uncovered  which  is 
significant  at  very  low  counting  rates.   In  our  system,  about  once  every  5  [min] ,  we 
witnessed  a  burst  of   50  to   100  counts.   This  means  about  a  1  percent  error  for  a 
counting  rate  of  30  [cs  ]  when  integrating  for  5  [min]   and  larger  errors  at  lower 
counting  rates.   This  effect  was  not  significant  down  to  290  [nm]   because  the  lowest 
counting  rate  at  this  and  longer  wavelengths  was  about  40  [cs   ].   The  effect  was 
significant  at  shorter  wavelengths  because  the  lowest  counting  rate  was  about  5  [cs  ] 
at  289  [nm] ,   2  at  288,  and  1  at  287.   However,  we  did  not  try  to  solve  this  problem 
because  scattering  was  a  serious  problem  below  290  [nm] ,   and,  as  you  will  see  shortly, 
we  were  not  able  to  correct  for  this  scattering.   So  an  additional  error  was  immaterial. 
We  were  only  able  to  satisfy  our  desired  overall   5-10  percent  uncertainty  at  wave- 
lengths of  290  [nm]   and  higher. 

Relative  to  spectral  distribution  —  wavelength  accuracy.   The  most  obvious  concern 
relative  to  the  large  difference  in  spectral  distribution  between  the  standard  source  and 
the  sun  is  that  of  wavelength  accuracy.   The  spectral  irradiance  of  the  sun  falls  off  at 
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a  rate  of  about  1  percent  per   0.01  [nm] .   Therefore  a  wavelength  uncertainty  of   0.01  [nm] 
or  less  is  desired. 

The  first  step  in  characterizing  the  spectroradiometer  relative  to  wavelength  was  to 
check  its  reproducibility.   The   296.728  [nm]   mercury  spectral  line  from  a  low  pressure 
mercury  source  was  used  for  this  purpose.   With  the  spectral  slit  width  of  the  instrument 
being   0.2  [nm] ,   the  "peak"  of  the  line  was  too  broad  to  utilize  its  position  to  check 
wavelength  reproducibility  to   0.01  [nm] .   A  wavelength  reference  which  is  orders  of  mag- 
nitude more  sensitive  than  the  peak  of  a  spectral  line  is  the  steep  sides  of  the  line. 
When  the  spectral  slit  width  is  large  compared  to  the  spectral  line  width,  this  is  really 

the  side  of  the  slit-scattering  function.   A  suitable  region  to  use  is  that  part  of  the 

_2 
spectrally  scanned  line  which  is  a  factor  of   10    less  than  the  peak.   The  rate  of  change 

_2 
of  signal  with  wavelength  at   10    of  the  peak  is  at  least   50  percent   per   0.01  [nm] 

compared  to  about   0.5  percent  per  0.01  [nm]   at  the  peak. 

We  found  that  the  wavelength  reproducibility  for  our  monochromator  varied  from  0.01 
[nm]   to  about   0.05  [nm]   depending  on  the  prior  (immediate)  history  of  the  wavelength 
drive  and  setting.   The   0.01  [nm]   figure  relates  to  scanning  repeatedly  over  a  spectral 
line,  covering  a  range  of  about   1  [nm]   and  taking  a  few  minutes  for  each  scan.   The   0.05 
[nm]   reproducibility  corresponds  to  scanning  repeatedly  over  the  line  but  covering  a  range 
of  about   60  [nm]   or   30  [min]   for  each  scan.   By  cleaning,  tightening,  adjusting,  and 
lubricating,  it  was  possible  to  improve  the   0.05  [nm]   only  slightly  to   0.04  [nm] .   It 
was  clear  that  to  achieve   0.01  [nm]   wavelength  uncertainty  would  require  a  data-taking 
schedule  (short-term,  small  wavelength  range  repetitive  scans)  that  was  not  practical  in 
most  cases.   Another  possiblity  for  realizing  a   0.01  [nm]   reproducibility  was  to  use  the 
solar  structure  shown  in  figure  4  to  perform  a  wavelength  self-calibration  in  each  experi- 
mental run. 

In  addition  to  the  short-term  reproducibility  just  discussed,  wavelength  accuracy  may 
also  depend  on  somewhat  longer  term  changes  (few  hours)  caused  by  ambient  temperature 
variations.   The  effect  of  instrument  temperature  on  wavelength  setting  was  checked  and 
found  to  be  only  about   0.01  [nm]   per   3  [°C]   change  in  temperature.   This  was  considered 
to  be  an  insignificant  factor  in  our  case  compared  to  the  short-term  reproducibility. 

Wavelength  accuracy  is  also  affected  by  any  lack  of  symmetry  in  the  slit-scattering 
function,  particularly  when  the  wavelength  determination  is  made  by  the  technique  (sides  of 
slit-scattering  function)  described  earlier.   This  will  be  discussed  in  the  next  section 
along  with  the  effects  of  distortion  and  scattering. 

Relative  to  spectral  distribution  —  distortion.    The  central  portion  of  the  slit-scat- 
tering function  is  required  in  order  to  ascertain  and  correct  for  the  amount  of  distortion. 

-2 
We  determined  the  central  portion  of  the  slit-scattering  function  down  to  about   10    of 

its  peak  value  by  spectrally  scanning  over  various  mercury  lines  between   289.4  and  365 

_2 
[nm] .    It  was  not  possible  to  get  lower  than  about   10    with  a  mercury  line  because  of 
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the  background  radiation  from  other  lines,  particularly  the  254  [nm]   line.   For  values 

-2 
below  10   ,   we  used  a  laser  line  at   290  [nm] .   We  obtained  this  UV  line  from  a  tunable 

CW  dye  laser  operating  at   580  [nm]   by  frequency  doubling  with  an  appropriate  crystal. 
The  laser  beam  was  positioned  along  the  optic  axis  of  the  monochromator  with  the  diffuser 
in  place.   The  size  of  the  beam  was  adjusted  so  that  the  instrument  was  irradiated  with 
about  the  same  solid  angle  as  that  subtended  by  the  sun.   A  spectral  scan  over  these  lines 
should  be  a  good  approximation  to  the  central  portion  of  the  slit-scattering  function  (p. 
20,  910-4)  [6].   When  irradiated  on  axis,  the  central  portion  was  quite  asymmetric.   How- 
ever, the  degree  of  asymmetry  changed  with  the  direction  of  the  beam  entering  the  mono- 
chromator.  It  was  possible  to  find  a  particular  direction  slightly  off  axis,  for  which 
good  symmetry  was  realized.   This  was  done  and  all  solar  measurements  made  without  the 
diffuser  were  made  with  the  solar  beam  incident  from  this  direction.   The  asymmetry  was 
also  quite  pronounced  in  the  diffuse  mode;  that  is,  with  the  diffuser  filled  with  optical 
radiation.   However,  by  trial  and  error,  it  was  possible  to  find  an  area  on  the  diffuser 
(equal  to  about  one-half  the  total  area)  which  also  resulted  in  a  symmetric  slit-scattering 
function. 

Figure  14  shows  an  example  at   296  [nm]   of  that  portion  of  the  slit-scattering  func- 
tion obtained  with  the  mercury  lines.   Also  shown  are  triangular  and  Gaussian  functions 
having  the  same  spectral  slit  width.   These  functions  do  a  fair  job  of  representing  our 
measured  slit-scattering  function  down  to  about  0.3  of  its  value.   From  0.3   to  about 
10   ,   a  much  better  representation  is  the  following  function: 


z(A  ,A)   =  exp[-126.38-(A  -A)2]  +h   exp[-2. 983-10  3-(A  -A)  2 


-1.2037-104-(Ao-A)4], 


(1.17) 


where  A   is  the  monochromator  setting  and   A   is  the  wavelength  in  nanometers  [nm] .   The 

°  -5 

slit-scattering  function  for  values  less  than  10  "   is  discussed  in  the  next  section. 

Incidentally,  the  spectral  slit  width  varied  slightly  with  wavelength  from  0.183  [nm]   at 

289  [nm]   to   0.165  [nm]   at   312  [nm]   and  then  back  to   0.176  [nm]   at   334  and  365  [nm] . 

Convolution  calculations  indicated  that  it  was  more  important,  relative  to  distortion 
effects,  to  know  that  the  slit-scattering  function  was  symmetric  than  to  know  its  exact 
shape  or  width.   Also,  the  existence  of  symmetry  meant  that  wavelength  calibrations  could 

be  performed  by  determining  the  average  of  the  wavelengths  on  the  two  sides  of  the  slit- 

-2 
scattering  function  where  the  value  of  the  function  is  about   10    of  its  peak  value. 

Relative  to  spectral  distribution  —  scattering.    In  order  to  correct  for  scattered  radia- 
tion that  we  anticipate  to  be  significant  below  about   292  [nm] ,   the  long  wavelength  wing 
of  the  slit-scattering  function  is  required.   We  tried  to  obtain  these  wings  with  our  tuned 
UV  laser,  but  the  signal-to-noise  ratio  was  too  small  to  get  lower  than  about  2  x  10   , 
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Figure  14.   Measured  slit-scattering  function  z(A  ,X)   at   X    296.62  [nm] 

and  triangular  and  Gaussian  functions  having  the  same  spectral  width. 


29 


even  with  an  integrating  time  of  one  hour.   The   2  x  10    was  the  value  at  about   1.8 
[nm]   from  the  peak.   We  obtained  a  rough  idea  of  the  magnitude  of  the  more  distant  wing 
by  scanning  over  the  short  wavelength  portion  (p.  22,  910-4)  [6]  of  the  253.7  [nm] 
spectral  line  from  a  very  intense  low  pressure  mercury  source  using  suitable  filters  to 

remove  radiation  other  than  that  scattered  from  the   253.7  [nm]   line  itself.   At   200 

—8 

[nm]   the  result  was   1.4  x  10   .    Thus  the  slit-scattering  function  at   +53.7  [nm] 

—8 
(i.e.,  253.7-200)   from  the  peak  is  about  10    when  the  monochromator  is  irradiated 

with  highly  incoherent  radiation  that  fills  the  f/5  field  of  view. 

So,  much  to  our  disappointment,  we  were  not  able  to  obtain  an  accurate   (5  or  10 
percent)   determination  of  the  distant  wing  of  the  slit-scattering  function  and  would 
not  be  able  to  accurately  correct  for  Scattering.   This  will  be  possible  when  a  tunable 
laser  is  available  with  about   10  [mW]   of  power  at   290  [nm] .    Until  then,  the  short 
wavelength  limit  for  our  accurate  measurements  will  be  determined  by  the  wavelength  at 
which  scattering  begins  to  be  significant.   The  solar  data  itself  will  reveal  this 
wavelength,  as  will  be  shown  shortly. 

DATA  TAKING.    In  this  section  we  discuss  not  only  details  of  the  actual  data  taking  but 
also  alignment  of  the  monochromator  relative  to  the  sun,  alignment  of  the  tungsten  lamp 
standard,  and  calibration  using  the  standard. 

Monochromator  alignment.    Careful  alignment  of  the  monochromator  was  essential  prior  to 
making  any  solar  measurements,  particularly  when  the  diff user  was  not  used.   Figure  8 
shows  that  an  angular  misalignment  of  a  degree  would  result  in  a  significant  error. 
Also,  as  mentioned  in  the  previous  section,  an  angular  misalignment  would  destroy  the 
spectral  symmetry  of  the  slit-scattering  function.,  reducing  wavelength  accuracy  and 
increasing  distortion. 

The  monochromator  was  aligned  in  the  following  manner.   A  low-power  He/Ne  laser  was 
set  up  in  the  laboratory  and  aligned  so  that  its  beam  passed  upward  approximately  through 
the  center  of  the  solar  tracker.   Then  the  laser  was  fine-adjusted  in  position  and 
orientation  so  that  its  beam  (see  figure  15)  did  not  move  relative  to  the  upper  aperture 
when  the  tracker  followed  the  sun.   This  meant  that  the  laser  beam  was  striking  each 
tracker  mirror  at  the  intersection  of  the  plane  of  the  mirror  with  the  vertical  and 
horizontal  axes  of  rotation,  as  indicated  in  figure  15.   After  this  condition  was 
realized,  the  upper  aperture  was  positioned  and  secured  so  that  the  laser  beam  passed 
through  its  center.   The   5  [in]   diameter  lower  aperture  was  also  positioned  so  that 
the  laser  beam  passed  through  its  center.   Next  the  laser  was  set  up  on  the  roof  and 
again  aligned  so  that  its  beam  went  downward  through  the  centers  of  both  apertures.   The 
monochromator  was  then  positioned  and  oriented  so  that  the  laser  beam  struck  the  center 
of  the  entrance  slit  and  coincided  with  the  optic  axis.   This  was  easy  to  accomplish 
because  the  tube  of  the  dif fuser/f ilter  system  had  been  made  and  mounted  so  that  its  top 
surface  was  perpendicular  to  the  optic  axis.   A  microscope  cover  glass  placed  on  this 
surface  was  therefore  perpendicular  to  the  optic  axis  and  the  monochromator  was  rotated 
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Figure  15.   Use  of  a  laser  beam  to  find  the  axes  of  the  solar  tracker. 
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until  the  laser  beam  retroref lected  back  on  itself  (or  the  center  of  the  5  [in] 
aperture).   Finally,  the  monochromator  was  rotated  about  its  entrance  slit  by  the  amount 
required   (0.4  [deg])   for  a  spectrally  symmetric  slit-scattering  function. 

In  performing  the  above  alignments  it  was  essential  to  use  a  heavy,  stable,  laser 
mount  with  micrometer  adjustments  for  translation  and  two  orthogonal  rotations  about 
axes  perpendicular  to  the  laser  beam.  In  addition,  it  was  convenient  to  use  circular 
aluminum  inserts  each  with  a  translucent  "button"  with  a  small  hole  drilled  through  the 
center  as  shown  in  figure  16.   The  insert  could  be  placed  in  the  aperture  and  used  as  a 
convenient  indicator  of  the  aperture  center  while  adjusting  the  laser  beam. 


CIRCULAR  ALUMINUM  INSERT 


TRANSLUCENT  BUTTON 
WITH  CENTER  HOLE 


W///////////M  Wmzmzzzzp 


uuunuuuuus 


szzznznzmzzm 


APERTURE 


Figure  16.    Inserts  used  to  indicate  the  centers  of  the  tracker  apertures. 
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Alignment  of  standard.    For  the  spectroradiometer  calibration  it  was  necessary  to  posi- 
tion the  standard  lamp  precisely.   The  procedure  was  as  follows.   The  He/Ne  laser  was 
mounted  on  the  roof  above  the  tracker  and  aligned  until  the  laser  beam  went  through  the 
center  of  the  button  in  the  upper  aperture  insert  (see  figure  16)  and  was  incident  on  the 
center  of  the  monochromator  entrance  slit.   Another  insert  similar  to  that  in  figure  16 
was  fitted  into  a  rectangular  piece  of  aluminum  and  positioned  in  the  lower  aperture  as 
indicated  in  figure  17  so  that  the  laser  beam  traversed  its  button  also.   The  laser  was 
then  set  up  in  the  laboratory  and  aligned  so  that  its  beam  traversed  the  buttons  in  the 
lower  and  upper  apertures  in  the  opposite  direction.   Then,  the  tungsten  standard  was 
mounted  on  the  roof  tracker  so  that  its  center  was  irradiated  by  the  laser  beam.   The  lamp 
mount  was  an  aluminum  tubular  structure  of  sufficient  lightness  pot  to  limit  the  smooth 
motions  of  the  tracker  and  of  sufficient  rigidity  so  that  the  lamp  was  stably  positioned 
relative  to  the  3.25  [in]   tracker  aperture.   The  distance  between  the  entrance  slit  and 
the  lamp  was  selected,  as  mentioned  earlier,  so  that  the  lamp's  largest  dimension  sub- 
tended about   0.5  [deg]   (same  as  the  sun)  at  the  entrance  slit.   The  actual  distance  was 
measured  to  be  396.3  [cm]   with  an  uncertainty  of  about  0.5  [cm]. 


TRANSLUCENT  BUTTON 


5   [in]  LOWER  APERTURE 


ALUMINUM  RECTANGLE 


ALUMINUM  INSERT 


Figure  17.   Looking  down  from  mirror  1  in  tracker  at  aluminum 
pieces  and  button  used  in  lamp  adjustment. 
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Spectral  irradiance  calibration.    The  standard  lamp  calibrations  had  to  be  frequent 
enough  and  extensive  enough  relative  to  wavelength,  tracker  position  and  filters  so  that  a 
sufficiently  accurate  calibration  was  available  for  any  time,  wavelength  and  measurement 
configuration  (tracker  position  or  filter)  for  which  solar  data  was  taken. 

Relative  to  time,  calibrations  were  performed  as  often  as  practical,  primarily  when- 
ever weather  permitted.   This  turned  out  to  be  at  least  every  other  night  during  the 
period  in  which  good  solar  data  were  obtained.   Successive  calibrations  taken  about 
48  hours  apart  differed  by  a  maximum  of  a  few  percent.   Relative  to  wavelength,  standard 
lamp  calibrations  were  obtained  between  340  and  280  [nm]   at  5  [nm]   intervals.   Previous 
experience  had  shown  this  to  result  in  a  least  squares  fit  with  a  standard  deviation  of 
the  fit  of  a  few  tenths  of  a  percent.   To  confirm  that   5  [nm]   intervals  were  adequate, 
one  run  was  also  taken  at   1  [nm]   intervals.   To  obtain  the  effect  of  tracker  position, 
the  340-280  [nm]   runs  were  obtained  for  about  every  half  hour  of  tracker  position  between 
10:00  and  16:30  EST.   Figure  9  is  a  typical  example  of  a  relative  tracker  calibration. 
Though  a  complete  calibration  relative  to  the  tracker  was  not  performed  every  evening 
during  which  calibrations  were  performed,  sufficient  data  were  taken  to  show  any  changes 
in  the  relative  tracker  calibration  over  the  eight  days  of  data  taking.   No  significant 
change  was  detected.   Sufficient  standard  lamp  data  were  also  taken  with  the  neutral 
density  and  solar  blind  filters  and  the  diffuser  to  determine  their  transmittances  and  to 
detect  any  changes  in  their  transmittances  over  the  eight  days  of  operation.   None  were 
observed.   The  spectral  transmittances  of  each  of  these  optical  elements  were  obtained  in 
the  wavelength  range  in  which  they  were  used  (neutral  density  filter  and  diffuser,   340- 
305  [nm] ;   solar  blind,   296-287  [nm]).    For  optimum  accuracy,  it  was  necessary  to 
obtain  this  transmittance  with  the  same  measurement  beam  configuration  that  was  used  in 
taking  the  solar  data  (as  opposed,  e.g.,  to  using  a  general  purpose  spectrophotometer  to 
determine  the  transmittances) .   This  eliminated  the  need  for  a  tedious  characterization  of 
the  filters  and  the  optical  beam  relative  to  uniformity  and  direction. 

Wavelength  calibration.    Wavelength  calibrations  were  made  every  morning  and  late  afternoon. 
Five  mercury  spectral  lines  at   289.360,  296.728,  312.566,  334.148,  and  365.015  [nm]   from 
a  low  pressure  pen  ray  lamp  were  used.   The  lamp  was  used  with  a  5  [mm]   aperture  in 
front  of  it  and  with  the  aperture  at  a  distance  of  about   573  [mm]   from  the  slit,  so  that 
the  largest  plane  angle  subtended  by  that  aperture  at  the  slit  was  0.5  [deg]  .   The  mercury 
lamp  and  its  aperture  were  also  positioned  so  that  the  path  of  this  beam  through  the 
monochromator  was  the  same  as  that  of  the  solar  beam.   As  discussed  earlier,  the  wave- 
length encoder  readout  for  each  scanned  spectral  line  was  taken  to  be  the  average  of  the 

-2 
readouts  on  each  side  of  the  peak  where  the  signal  was  exactly  equal  to   10    that  of  the 

peak   (10    for  the  289.36  [nm]   line).   To  obtain  an  estimate  of  the  wavelength  hysteresis 

effects  that  had  been  observed,  wavelength  calibrations  were  performed  starting  at  about 

289  [nm]   and  going  to  about   365.4  [nm]   and  then  back  to   289  [nm] . 

Wavelengths  for  each  of  the  five  lines  were  always  obtained  in  a  spectral  scan  that 
was  decreasing  in  wavelength.   This  was  the  case  even  when  performing  calibrations  starting 
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with  the   289  [nm]   line  and  ending  with  the   365  [nm]   line.   Each  line  in  this  sequence 
was  still  scanned  in  decreasing  wavelength  to  minimize  backlash  errors.   A  typical  cali- 
bration correction  curve  is  shown  in  figure  18  where  the  difference  between  the  published 
wavelength  of  a  line  and  the  corresponding  wavelength  encoder  readout  is  plotted  against 
wavelength.   The  largest  wavelength  hysteresis  usually  observed  was   0.04  [nm] .    The 
dashed  line  in  the  figure,  the  mean  of  the  calibrations  in  the  two  wavelength  directions 
(indicated  by  arrows),  was  used  as  the  calibration  for  the  time  these  wavelength  data  were 
obtained. 

Solar  measurements.    Solar  data  were  taken  in  Gainesville,  Florida  between  June  11,  1980 
and  June  18,  1980  whenever  the  sky  was  sufficiently  clear.   This  meant  that  there  were  no 
visible  clouds  in  the  sun's  line  of  sight.   Four  types  of  data  were  obtained  for  four 
different  purposes.   These  were   (1)  short  wavelength  cut-off  data   (287-296  [nm]   at 
1  [nm]   intervals),   (2)  data  for  ozone  determinations   (294.7  -  305.2  [nm]   at  a  suffi- 
cient number  of  wavelengths  so  that  the  spectral  structure  could  be  used  to  obtain  at 
least  5  spectral  irradiances  for  which  the  wavelength  uncertainty  was  not  greater  than 
0.01  [nm]),   (3)  data  for  ozone  determination  at  the  "Dobson  wavelengths"   (305.5  -  340) 
[nm]   and   (4)  data  at  various  wavelength  intervals  of   0.1  [nm]   to   2  [nm]   over  the 
wider  spectral  interval  of  the  solar  spectral  irradiance  which  would  be  required  for 
deconvolution  calculations   (295  -  340)  [nm] .    In  all,  eighty-one  different  experimental 
runs  were  carried  out. 

The  electronic  components  of  the  spectroradiometer  were  left  on  continuously  for  the 
entire  8  days  of  data  taking.   The  anode  to  cathode  voltage  on  the  photomultiplier  was 
maintained  at   1500  [V] .    The  discriminator  voltage  was  set  to  produce  a  dark  count  rate 
of  about   0.3  [c«s   ]  . 

Before  starting  an  experimental  run,  a  check  was  made  to  insure  that  the  shadow  of 
the  3.25  [in]   tracker  aperture  was  concentric  with  the  dif fuser/filter  system.   If  the 
shadow  was  not  centered,  it  was  made  so  by  overriding  the  computer  control  of  the  tracker 
stepping  motors  with  a  manually  operated  control.   Next,  the  appropriate  filter  or  dif- 
fuser  cap  was  inserted  and  a  computer  program  selected  for  one  of  the  four  types  of  data 
runs.   In  order  to  observe,  during  each  run,  all  parts  of  the  spectrum  under  closely 
similar  atmospheric  conditions  and  solar  paths,  integration  times  were  limited  to  about 
1  [s]   or  to  that  required  to  obtain  about   10,000  [c] .   An  experimental  run  required 
between   3  and  30  [min]   with  most  of  them  about   5  [min] .    In  order  to  minimize  this  run 
time,  data  being  taken  were  usually  stored  on  the  computer  disk,  and  a  hard  copy  made  at 
end  of  day.   Continuous  visual  monitoring  of  the  sun  was  maintained  during  the  run  to 
insure  that  no  visible  clouds  were  in  line  of  sight  of  the  sun.   Notes  were  taken  during 
the  experimental  run  describing  environmental  conditions  and  commenting  on  anything  that 
occurred  that  was  considered  potentially  useful.   As  many  of  the  measurements  as  possible 
were  made  within  one  hour  of  high  noon  (^12:30  EST).   When  comparisons  were  to  be  made; 
e.g.,  comparing  the  Dobson  and  "Sharp  Cut-Off"  ozone  determinations,  one  type  of  run  was 
made  twice,  and  it  bracketed  the  other  type  in  time;  and  all  runs  were  made  in  as  short  a 
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Figure  18.    Typical  correction  curve  for  the  monochromator 

wavelength  readout.   Maximum  observed  wavelength 
hysteresis  was  about  0.04  [nm] . 
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time  as  possible.   In  this  way,  atmospheric  changes  were  less  likely  to  occur  during  the 
comparison,  and  if  they  did  occur,  would  be  more  likely  to  be  detected.   Data  taking  was 
automatically  interrupted  via  the  computer  program  when  a  filter  or  diffuser  had  to  be 
inserted.   Ambient  temperature  control  was  limited  to  manual  control  of  window  air  condi- 
tioners but  the  air  temperature  near  the  monochromator  was  taken  and  automatically  recorded 
about  once  per  minute.   Barometric  pressure  readings  were  obtained  continuously.   Records 
for  the  weather  parameters  such  as  humidity,  visibility,  temperature,  and  wind,  taken  each 
hour,  were  available  from  the  Gainesville  Weather  Station.   Table  1  is  an  example  of  a 
typical  computer  printout  of  data  obtained  in  one  short  run. 


Table  1.  Typical  computer  printout  for  a  short  run  of  data. 


DATA  FILE  NAME-  BrWSOLBl 
DATE  6/18/80 

TEMPERATURE  2  5.03  9  [°C] 


TIME  11:51:25  EST 


WAVELENGTH 

296.002 
295.000 
294.001 
293.002 
292.002 

TIME-  11:53:52  EST 

291.000 
289.998 
288.998 
287.999 
287.000 

TIME-    12:09:25    EST 


PHOTON    COUNT 

319828 

93692 

52306 
256625 
119052 


COUNTING    RATE       [cps] 

159914.00 

46846.00 

26153.00 

8554.17 

3968.40 


TEMPERATURE    25.097  [°C] 

26241  874.70 

17264  287.73 

6177  51.48 

1278  10.65 

3029  5.05 

TEMPERATURE-    23.030     [°C] 


DATA  REDUCTION  and  UNCERTAINTIES.        Data  reduction  means  obtaining  the  solar   spectral 

irradiance     E,      at   the  wavelengths   of   interest   or   the  average  of  E       weighted  by   the   slit- 
A  A 

scattering   function.      In   this  problem,    because   the   spectral   slit  width  is   so  much  greater 
(0.18    [nm])      than  the  wavelength  differences   for  which  large   changes   occur   in   the   solar 
spectrum      (<0.01    [nm]),      we  must   be   satisfied  with  a  weighted   average  of     E    . 
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In  reducing  and  analyzing  the  data,  we  start  with  the  measurement  equation,  eq. 
(1.12), 


S(AQ,t,T)       C(S').S'Uo,t,T)         J   Ex(X,t).RE(Xo,A,t,T).d; 

_  j  AA 

|   E^(A,t).RE(Ao,A,t,Ts).d; 

J  A) 


'AA  A       _  " (1.18) 


SS(AQ,t,Ts)       C(S,S)-S'S(Ao,t,Ts) 


where  the  independent  variables   A    (wavelength  setting  of  the  monochromator) ,  t   (time 
at  which  S'   was  obtained),  and  T   (tracker  orientation  indicated  by  time)  have  been 
inserted  to  emphasize  that  the  responsivity  is  a  function  of  these  variables.   The  irradi- 
ance  responsivity  may  be  written 

RE(Ao,A,t,T)   =  T(A,t)-z(Ao,A,T).R^(Ao,t,T),  (1.19) 

where  t(A,t)   is  the  spectral  transmittance  at  time   t  of  any  filter,  including  the 

diffuser,   z(A  ,A,T)   is  the  slit-scattering  function  (p.  15,  910-4)  [6]1,  and  R_(A  ,t,T) 
o  bo 

is  the  responsivity  factor  for  tracker  orientation  T.   When  we  substitute  eq.  (1.19)  into 
the  measurement  equation,  eq.  (1.18),  we  have 

S'(Ao,t,T)       C(S'S).R^(Ao,t,T).|   EA(A,t).x(A,t).z(Ao,A,T)-dA 


AA 


S,S(A   t,T  )      C(S')-r£(A  ,t,T  ).[   E,(A,t).z(A  ,A,T  ) -d) 
os  E  o    s    . ,  A         os 

1  AA 


(1.20) 


g 

where  x   for  the  case  of  the  calibration  signal,   S'  ,   has  been  set  equal  to  one  because 
no  filters  were  normally  used  in  obtaining   S'  .   The  unknowns  in  this  equation  are 

E. (A,t)   and  z.   We  have  reliable  measurements  for  the  central  part  of   z,   approximately 

A 

±0.4  [nm]   on  either  side  of  the  peak,  and  we  have  an  upper  limit  to  the  magnitude  of  z 

—8 
at  +50  [nm]   from  the  peak.   This  is  about   10   .    Roughly  sketching  in  the  intermediate 

values  of  z  and  then  using  this  slit-scattering  function  in  a  simple  deconvolution  (p. 

38,  910-4)  [11]  on  spectral  runs  for  which  data  were  obtained  to   287  [nm] ,   we  found  that 

the  computed  scattering  was  about  100  times  more  than  the  observed  signal  count  rate   S' 

at  287  [nm] .   The  only  explanation  we  have  for  this  huge  discrepancy  is  that  the  long 


*In  TN  910-4  [6],  the  slit-scattering  function  has  the  functional  form  z(A  -A).    It  is 
not  necessary  to  restrict  it  to  that  form  in  this  investigation,  and  we  define  it  here 
more  generally  as 


RE(A,AQ,t,T) 
R|(Xo,t,T) 


z(A  ,A,T)   = 


38 


wavelength  wing  of  the  slit-scattering  function  for  our  instrument  when  observing  the  sun 

—  1 0  —8 

without  a  diffuser1  is  about   10     rather  than   10   .   Calculations  show  that  with  a 

slit-scattering  function  equal  to   10     in  the  long  wavelength  wing,  scattering  would  be 
negligible  down  to   290  [nm] .    This  is  experimentally  confirmed  by  comparing   S'(287), 
which  is  about  equal  to  or  greater  than  the  scattered  signal  at  287  or  290  [nm] ,   to 
S'(290).   Since  there  is  not  enough  power  in  the  UV  in  available  CW  lasers  to  determine 
the  distant  wings  of  z   accurately,  the  lower  wavelength  limit  for  accurate  data  reduc- 
tion is  about   290  [nm] .    And  since  scattering  is  insignificant  for   X  >  290  [nm]  ,  we 
may  limit  AX   to  only  the  central  portion  of  z.    Trial  convolutions  showed  that  it  was 
sufficient  for  AX   to  be  equal  to   0.44  [nm] ;   that  is,   0.22  [nm]   on  either  side  of  the 
peak  of  z. 

For  the  wavelength  range   X  ±0.22  [nm] ,   eq.  (1.20)  can  be  simplified  as  follows 


X  +0.22 
o 


S'(X    ,t,T)  C(S,S).t(X    ,t)-R!(X    ,t,T)-  f    E.(X,t)-z(X    ,X,T)-d 

O  O  H       O  /  A  o 


X   -0.22 
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X  +0.22 
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X   -0.22  C(S')-S'(X    ,t,T).E^(X    ,t) 
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Xo+U'ZZ  A    ,°  C(S,S).S,S(X    ,t,T   )-t(X    ,t).P 

OS  o 


z(X    ,X,T   )-dX 
o  s 


X   -0.22 
o 

where  f 

R|(X    ,t,T.) 
p     =     _JL_° (1.23) 

is  the  polarization  correction  required  when  the  calibration  data  S'   are  obtained  with 

tracker  orientation  T   and  solar  data  S'   with  tracker  position  T.    Calculations 

s 

involving  eqs.  (1.20)  and  (1.21)  show  that  this  simplification  produces  a  negligible  error 
(<0.3  percent).    From  previous  experience,  we  expect  that  the  difference  between  z(X  ,X,T 


because  the  terrestrial  sun  is  confined  to  such  a  small  angle  (0.5  [deg]),  its  irradiance 
on  the  entrance  slit  of  the  spectroradiometer  is  highly  coherent.   Calculations  of  the 
type  described  by  Shumaker  [26]  show  that  for  the  f/5  spectroradiometer  that  we  used  the 
distant  wings  of  the  slit-scattering  function  could  be  reduced  by  this  order  of  magnitude. 
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and  z(A  ,  A,T  )   for  A  ±0.22  [nm]   is  small,  and  we  assume  the  effect  of  this  difference 
in  eq.  (1.22)  is  insignificant. 

All  the  data  and  the  equations  resulting  from  the  various  fits  were  stored  on  floppy 

disks.   It  was  therefore  a  simple  matter  to  write  a  program  to  compute  E,  (A  ,t).   This 

-z  X  ° 

was  done  and  E   was  computed  for  all  wavelengths  in  the  81  experimental  runs.   Examples 

A 

of  the  data  and  results  printed  out  from  these  computations  are  shown  in  tables  B-l  through 
B-7,  in  Appendix  IB.   Keep  in  mind  that  these  results  have  not  been  corrected  for  scat- 
tering.  However,  we  believe  this  to  be  significant  only  below  290  [nm] .   Also,  they 
have  not  been  corrected  for  possible  differences  between  using  and  not  using  a  diffuser, 
which  will  be  addressed  later. 

In  order  to  estimate  the  uncertainty  in  our  determination  of  E    [eq.  (1.22)],  we 

A 

shall  address  the  way  in  which  each  independent  and  dependent  variable  in  this  simplified 
equation  was  obtained. 

The  first  variable  is   A  ,   which  is  the  (corrected)  wavelength  setting  of  the  spectro — 
radiometer.   This  is  actually  the  peak  or  center  wavelength  of  our  slit-scattering  function. 
The  wavelength  calibrations  that  were  taken  every  morning  and  late  afternoon,  one  of  which 
is  shown  in  figure  18,  were  generally  used  to  obtain   A  .   The  data  represented  by  the 
dashed  curve  in  the  figure  was  fitted  to  a  cubic  equation  in  A .    A  linear  interpolation 
in  time  between  the  morning  and  afternoon  calibrations  was  used  to  determine  the  wave- 
length correction  at  any  time  t.    These  wavelength  calibrations  were  estimated  to  have 
an  uncertainty  ranging  from  ±0.02  [nm]  at  290  [nm] ,   ±0.03  [nm]  at  295  [nm] ,   ±0.03  [nm] 
at  305  [nm]   to  ±0.02  [nm]  at  340  [nm] .    A  smaller  wavelength  uncertainty,  estimated  to 
be  ±0.01  [nm] ,   was  obtained  in  the  data  taken  for  ozone  determinations  by  using  the 
spectral  structure  shown  in  figure  4.   In  this  case,  during  a  solar  run,  about  five  data 
points  were  taken  at   0.02  [nm]   intervals  through  5  or  6  of  the  troughs  in  figure  4.   The 

observed  trough  wavelengths  were  then  compared  with  those  obtained  by  convolution  of 

— z 
E  (A,t),   where  E  (A,t)   was  obtained  in  the  manner  described  in  Appendix  A.   The  E 

A  A  A 

uncertainties  associated  with  these  various  wavelength  uncertainties  vary  with  wavelength 
and  are  shown  in  table  2  along  with  the  other  uncertainties  to  be  discussed  in  this  section. 

.  The  time  (Eastern  Standard),   t,   was  provided  by  a  clock  in  the  computer  and  was 
accurate  to  about  one  second.   Time  was  stored  along  with  the  data  at  intervals  of  about 
one  minute.   The  time  parameter  was  used  to  compute  the  orientation  of  the  tracker  for 
polarization  corrections  and  for  linear  interpolations  of  the  standard  lamp  calibration, 

S'    and  introduced  negligible  uncertainty  in  E..  . 

A 

The  signal  counting  rate,   S',   for  solar  photons  was  obtained  with  an  integration 
time  varying  from  one  second  to  10  minutes,  chosen  to  realize  a  standard  deviation  for  S' 
of   1  percent  or  less.   In  photon  counting,  the  standard  deviation  of  the  signal  is 
approximately  equal  to  the  square  root  of  the  total  counts.   The  dark  count  rate  was 
always  set  to  about   0.3  [cs   ]   by  adjusting  the  discriminator  voltage  and  was  sub- 
tracted from  the  observed  signal  rate  when  significant.   Time  intervals  for  integrations 
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were  determined  by  the  computer  with  negligible  uncertainty.   Time  was  the  basis  for 
automatic  control  of  the  tracker  orientation  such  that  the  solar  beam  remained  centered 
relative  to  the  entrance  slit  of  the  spectroradiotneter.  This  automatic  mechanism  mal- 
functioned rather  frequently  during  the  runs,  requiring  manual  control.  We  have  assigned 
a  2  percent  uncertainty  for  this  malfunctioning  as  indicated  in  table  2. 

The  correction  C(S')   (for  not  being  able  to  resolve  all  the  photon  counts)  given  in 
eq.  (1.15b)  is  considered  uncertain  by  about  1  percent  for  counting  rates  from 
5  x  10  [cs  ]   to  the  highest  observed.  At  lower  counting  rates,  its  uncertainty  is 
estimated  to  decrease  monotonically  to  a  few  tenths  of  one  percent  at  10   [cs  ]. 

The  spectral  irradiance  of  the  1000  [W]  frosted  tungsten  halogen  standard  Ef(X,t) 

A 

was  obtained  by  a  comparison  calibration  with  three  NBS  spectral  irradiance  1000  [W] 
standards  before  and  after  the  eight-day  period  of  solar  observations.   This  was  performed 
at  the  normal  distance  of  50  [cm]   and  at  5  [nm]   intervals  for  the  wavelength  range 
280-350  [nm] .   No  significant  change  was  observed  in  the  two  calibrations  obtained  8  days 
apart.   The  spectral  irradiance  of  the  frosted  lamp  at  the  distance  used,   396.3  [cm] 
without  diffuser  and  384.8  [cm]  with  diffuser,  was  calculated  assuming  an  inverse 
square  law  dependence  relative  to  50  [cm].   From  previous  experience,  the  uncertainty  of 
this  procedure  for  obtaining  E   of  this  type  lamp  at  the  above  distances  does  not  exceed 
0.5  percent.   In  order  to  obtain  the  spectral  irradiance  at  wavelengths  other  than  the 
5  [nm]   intervals  at  which  data  was  taken,  a  least  squares  fit  was  performed  to  the  expres- 
sion 

E*(X)   =   <AX  +  A2-X  +  A3-X2  +  A4-X3  +  A5.xVx"5-e[a+(bA)]  (1.24) 

s  2 

with  weights  of  1/(E?)   (NBS  Tech  Note  594-13)  [13].   If  only  a  linear  least  squares 

computer  program  is  available,  one  can  determine  the  quantities  a  and  b  from  a  fit  of 

Hn(E®.X5)  -  a  +  (b/X)  (1.25) 

with  weights  of  1.  Then  with  a  and  b  determined,  one  obtains  the  A   from  a  fit  of 

s  2  s 

eq.  (1.24)  with  weights  1/(E.).   The  uncertainty  of  E   at  the  distances  and  wave- 
lengths used  in  this  experiment  is  estimated  to  be  3  percent. 

g 
The  standard  lamp  calibration  data  S   were  obtained  as  a  function  of  wavelength  by 

fitting  an  equation  similar  to  eq.  (1.24)  but  without  the  cubic  and  quartic  terms  (these 

s  s 

terms  were  not  significant  in  the  S   fits)  and  with  a  weighting  of  l/(At«S  )  where  At 

is  the  integration  time  in  seconds.   In  linear  least  squares  fitting,  the  weights  are 

4 
proportional  to  the  reciprocals  of  the  variances.  Total  counts  varied  from  about  10  [c] 

at  280  [nm]   to  3  x  10   [c]   at  350  [nm]  when  an  integrating  time  of  5  [s]  was  used 

(no  filters  or  diffuser) .   Fits  performed  with  both  1  [nm]   and  5  [nm]   interval  data 
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demonstrated  that  the   5  [nm]   interval  spacing  was  sufficiently  close.   Successive  cali- 
brations taken  48  hours  apart  differed  by  a  maximum  of  about   2  percent  which  was  attrib- 
uted to  imprecision  (noise),  alignment  errors,  and  changes  in  the  tracker  mirrors.   The 
calibration  used  for  any  particular  solar  measurement  was  obtained  by  linear  interpolation, 
with  respect  to  time,  between  the  48-hour  calibrations.    The  uncertainty  of   S   was 
estimated  to  be   2  percent   as  indicated  in  table  2. 

Spectral  transmittance  data  of  the  solar  blind  filter  and  diffuser  were  fitted  with  a 
quadratic  function  of  wavelength  with  an  insignificant   (<0.5  percent)   standard  deviation 
for  the  fits.   The  solar  blind  filter  had  significant  non-uniformities  and  the  reproduci- 
bility in  positioning  this  filter  was  about  2  percent.   This  problem  did  not  exist  for  the 
diffuser.   The  neutral  density  filter  had  a  constant  transmittance  as  a  function  of  wave- 
length to  within  the  precision  of  the  measurements  of   1  percent.   It  did  not  possess 
significant  non-uniformities.   None  of  the  filters  showed  any  change  in  transmittance  over 
the  five  nights  of  calibration.   There  was  no  reason  to  expect  the  solar  blind  or  neutral 
density  filter  to  be  dependent  on  tracker  orientation,  and  this  was  confirmed  experimen- 
tally.  Since  the  output  signal  of  the  spectroradiometer  was  not  a  function  of  tracker 
orientation  with  the  diffuser  in  place,  but  without  it,  the  output  signal  did  vary  with 
tracker  orientation,  the  effective  transmittance  of  the  diffuser  was  always  determined 
relative  to  one  particular  tracker  orientation.   This  was  taken  to  be  the  orientation  at 
12:30  EST.   The  uncertainties  assigned  these  filters  were   3  percent   for  the  solar  blind 
filter,   2  percent   for  the  diffuser,  and   1  percent   for  the  neutral  density  filter. 

The  polarization  correction  for  tracker  orientation  was  obtained  by  two  quadratic 
spline1  "fits."   The  first  fit  was  used  to  interpolate  relative  to  tracker  orientation  for 
each  of  the   5  [nm]   lamp  calibration  data.   Since   Rp(^  >t,T  )   is  proportional  to 


*A  spline  "fit"  is  a  manner  of  interpolating  non-linearly  between  a  large  number  of  data 
points.   The  procedure  is  to  determine  the  coefficients  of  an  nth  order  (usually  n  =  2 
or  3)   polynomial  that  "goes  through"  the  first   n  +  1   data  points,  then  the  coefficients 
of  the  same  order  polynomial  that  "goes  through"  the  second  to   n  +  2   data  points,  and  so 
on  until  all  the  data  points  are  used.   Then  these  polynomials  are  used  to  calculate  the 
interpolated  values.   Since  two  or  more  polynomials  can  be  used  for  calculating  the  same 
interpolated  point,  a  standard  format  is  set  up  which  selects  a  particular  polynomial. 
For  example,  if   n  =  3,   a  typical  format  would  be  to  use  the  first  polynomial  for  inter- 
polating between  the  first   3   data  points  and  the  last  polynomial  for  interpolating 
between  the  last   3   data  points.   Each  of  the  remaining  polynomials  is  used  for  inter- 
polating between  the   2nd   and   3rd   point  used  in  determining  that  polynomial.   Once  a 
computer  program  has  been  written  to  accomplish  all  this,  it  is  a  simple,  fast  procedure 
to  carry  out.   We  have  often  used  spline  "fits"  when  we  wanted  data  at  equal  intervals  of 
the  independent  variable  but  the  data  wasn't  taken  at  these  intervals,  and  we  did  not  know 
an  appropriate  equation  to  use  for  a  least-squares  fit. 

43 


S'  (A  ,t,T  )  •  C(S'  ),   these  data  could  be  used  to  interpolate  relative  to  the  10:00  to 

OS 

17:00  tracker  orientations   (T) .   Curves  similar  to  that  shown  in  figure  9  were  computed 
and  stored  on  a  floppy  disk  for  each  of  the  5  [nm]   runs  from  340  to  280  [nm] .    The 
tracker  polarization  correction  P  eq.  (1.23)  was  determined  by  obtaining  the  ratio 

S,S(Xo,t,T) 


S,S(Ao,t,T') 


where  T  is  the  tracker  time  for  which  the  particular  solar  data  S'  was  obtained  and 
T'   is  the  closest  tracker  time  to  T  for  which  a  lamp  calibration  was  obtained.   These 
ratios  were  determined  for  each  of  the  5  [nm]   data  and  a  2nd  quadratic  spline  fit  was 
performed  as  a  function  of  wavelength,   The  uncertainty  of  P  was  estimated  to  be  2 
percent. 

In  order  to  detect  any  changes  or  drifts  with  respect  to  time  in  the  spectroradiometer 
responsivity  function,  the  standard  lamp  was  set  up  in  the  laboratory  and  observed  each 
morning  before  solar  measurements  began  and  in  the  late  afternoon  after  solar  measurements 
for  that  day  had  ended.   The  major  differences  between  these  lamp  observations  and  the 
roof  calibrations  were  that  the  former  did  not  include  the  tracker  and  did  not  have  to  be 
made  at  night  in  an  outdoor  environment.   These  made  the  "lamp  observations"  more  precise 
and  they  could  be  obtained  more  frequently.   In  order  that  corrections  C(S')   for  high 
count  rates  be  negligible,  the  neutral  density  filter  was  employed.   Figure  19  is  a  plot 
of  these  results  and  the  roof  calibrations   (S'  )   for   330,  335  and  340  [nm] .   The 
"inside-laboratory"  calibration  had  the  same  long-term  drift  as  the  roof  calibration. 
However,  short-term  changes  were  different.   We  do  not  know  the  reason  for  this,  but  it  is 
conceivably  due  to  laboratory  temperature  changes.   We  attributed  an  uncertainty  of   2 
percent   to  this  effect,  referring  to  it  in  table  2  as  "linear  interpolation  relative  to 
time  of  roof  calibration". 

As  mentioned  earlier,  we  had  planned  to  take  several  sets  of  data  at  the  longer 
wavelengths  with  and  without  a  diffuser  to  determine  directly  any  differences  between 
these  two  approaches.   The  set  of  runs  42,  43,  and  44  in  tables  B-5,  B-6,  and  B-7, 
respectively,  is  one  such  example.   All  these  comparisons  showed  that  there  was  a  difference 
of  about   5  percent,   independent  of  wavelength,  between  the  "diffuser"  and  "no-dif fuser" 
runs  with  the  latter  resulting  in  higher  values  of  spectral  irradiance.   We  have  not  been 
able  to  determine  the  cause  of  this  difference,  but  our  suspicion  is  that  it  is  due  to 
differences  in  the  beam  geometry  of  the  sun  and  standard  lamp  and  lack  of  sufficient 
uniformity  in  the  flux  responsivity.   We  believe  the  diffuse  data  is  correct  and  therefore 
must  reduce  all  the  data  taken  without  a  diffuser  (virtually  all  below  305  [nm])   by  the 
5.3  percent  difference  observed.   This  correction  has  not  yet  been  applied  to  the  non- 
dif fuser  data  in  tables  B-l  to  B-7  in  Appendix  IB.   The  uncertainty  of  such  a  correction 
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Figure  19.   Change  with  time  of  direct  mode  lamp  calibrations  (S'  ) 
for  300,   335,   and   345  [nm]   data. 
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has  been  estimated,  from  the  precision  of  the  dif fuser/no-dif fuser  comparisons,  to  be 
about   2  percent.   This  is  the  last  item  in  the  uncertainties  of  table  2. 

It  is  common  practice  to  make  an  estimate  of  the  overall  or  total  uncertainty  of 
experimental  measurements  by  combining  in  quadrature  (square  root  of  the  sum  of  the  squares) 
all  the  individual,  independent  uncertainties  such  as  those  listed  in  table  2.   This  has 
been  done,  and  the  result,  given  at  the  bottom  of  table  2,  varies  from   7.4  percent   at 
short  wavelengths  to   5.7  percent   at  long  wavelengths. 

There  are  applications  of  solar  terrestrial  spectral  irradiance  measurements  where 
only  the  ratio  of  values  at  wavelength  intervals  of   2  or  3  [nm]   are  needed.   Many  of  the 
uncertainties  in  table  2  are  significantly  less  for  such  relative  spectral  measurements, 
particularly  when  the  data  from  which  the  ratio  is  to  be  determined  have  been  obtained  in 
time  intervals  of  a  few  minutes.   We  estimate  that  the  total  uncertainty  for  such  relative 
spectral  measurements  is  about   3  percent   or  roughly  half  that  shown  in  table  2.   The 
component  uncertainties  contributing  to  this  smaller  uncertainty  are  wavelength — 2  per- 
cent,  noise  (S')  —  1  percent,   C(S') — 1  percent,   noise  (S'  ) — 1  percent,   E.. — 1  percent, 

A 

and  polarization — 1  percent. 

The  repeatability  of  non-diffuser  experimental  runs  taken  at  10  to  30  minute  intervals 
is  not  as  good  as  predicted  by  our  uncertainty  estimates.   Comparing  such  data  when  viewing 
conditions  were  good  and  wavelength  uncertainties  negligible,  we  found  systematic  dif- 
ferences that  averaged  4  percent   in  successive  runs.   We  did  not  see  this  large  a  dif- 
ference when  the  dif fuser  was  used.   However,  we  only  had  one  appropriate  set  of  dif fuser 
runs  for  checking  this  repeatability.   In  this  case,  there  was  a  difference  of  only 

1  percent.    The  only  factor  of  which  we  are  aware  that  could  cause  the  lack  of  repeata- 
bility of  4  percent   for  the  non-diffuser  runs  is  improper  alignment  relative  to  the  peak 
of  the  responsivity  function  (see  figure  8).   In  this  case,  systematic  errors  of  about 

1/2  deg   in  solar  tracking  during  a  run  would  account  for  the  results  observed.   The 
tracking  malfunction  uncertainty  in  table  2  addresses  this  type  of  behavior.   However,  the 
amount  assigned  this  uncertainty  (in  table  2)  accounts  for  only  half  the  non-repeatability 
observed,  and  we  wbuld  expect  this  uncertainty  to  reflect  itself  in  a  somewhat  less 
systematic  manner.   As  a  result  of  this  unaccountable  behavior  of  our  data,  we  add  another 

2  percent  uncertainty  to  the  totals  in  table  2  resulting  in  about   10  percent  at  the 
short  wavelengths  and  8  percent   at  the  long  wavelengths. 

The  question  arises  of  whether  or  not  a  significant  reduction  of  uncertainty  in  these 
solar  terrestrial  spectral  irradiance  measurements  is  likely  without  a  major  improvement 
in  the  spectroradiometer  and  standard.   We  believe  the  answer  is  yes  if  the  existing 
mirror  tracker  were  eliminated  and  the  monochromator  were  located  on  the  roof  and  itself 
mounted  on  a  very  accurate  solar  tracker.   Arrangements  also  would  have  to  be  made  for 
proper  temperature  control  of  the  spectroradiometer.   With  such  an  arrangement,  the  solar 
tracking  malfunction  uncertainty,  the  systematic  differences  referred  to  in  the  last 
paragraph,  the  polarization  correction  and  its  uncertainty,  and  half  the  uncertainty 
associated  with  the  standard  lamp  signal   S'    would  be  eliminated.   In  addition,  we 
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believe  the  uncertainty  associated  with  the  lack  of  uniformity  of  the  filters  and  with 
their  transmittance  measurements  could  be  reduced  by  a  factor  of  about  two.   Furthermore, 
whenever  possible,  one  should  use  a  diffuser.  This  is  no  problem  for  wavelengths  greater 
than  305  [nm] .   Between  295  and  305  [nm]   it  is  possible  only  when  integration  times  of 
100  to  10  [s],   respectively,  can  be  employed.   In  all,  it  should  be  possible  to  realize  a 
total  uncertainty  as  computed  in  table  2  of  about  6  percent  at  290  [nm]   and  4  percent 
between  about  300  and  340  [nm] . 

FUTURE  INVESTIGATIONS.   There  are  three  investigations  or  areas  of  research  to  which  we 
would  give  a  high  priority  in  any  future  efforts  in  UV  solar  terrestrial  radiometry. 

The  first  investigation  is  a  direct  test  of  the  effect  of  coherence  of  the  solar  flux 
on  spectral  scattering.  We  found  negligible  scattered  signal  at  290  [nm]  which  we 

attributed  to  the  high  coherence  of  the  flux  when  the  sun  itself  and  the  entrance  slit 

-4 
defined  the  solar  beam  (solid  angle  at  entrance  slit  of  approximately  10   [sr]).   This 

reason  for  the  low  scattering  could  be  confirmed  easily  by  first  spectrally  scanning 

between  295  and  290  [nm]  with  a  low  scattering,  double  monochromator  under  our  viewing 

conditions  and,  immediately  afterward,  performing  another  scan  with  an  appropriate  lens 

inserted  so  that  the  entrance  slit  was  irradiated  with  a  much  larger  solid  angle  (e.g., 

the  collimator  of  an  f/5  instrument  being  filled  or  a  solid  angle  of  approximately  0.03 

[sr]).   If  coherence  is  responsible  for  the  significantly  lower  scattering,  S(295)/S(290) 

from  the  two  scans  should  differ  by  a  factor  of  between  about  50  and  100. 

The  second  area  of  research  in  UV  solar  terrestrial  radiometry  that  we  believe  should 
receive  high  priority  is  extending  the  measurements  of  the  sun  itself  described  in  this 
chapter  to  measurements  of  the  sun  and  sky;  that  is,  hemispherical  viewing.   In  the  UV, 
the  solar  terrestrial  hemispherical  spectral  irradiance  is  about  twice  that  of  the  spec- 
tral irradiance  of  the  sun  itself  [14].  Ultraviolet  measurements  of  the  sun  itself,  as 
described  in  this  chapter,  are  of  major  interest  for  ozone  and  aerosol  determinations. 
Hemispherical  measurements  are  of  interest  in  studies  of  biological  effects.   The  usual 
way  in  which  hemispherical  measurements  are  made  is  with  a  diffuser  that  obeys  the  cosine 

law.  Development  and  accurate  characterization  of  a  suitable  diffuser  would  constitute 

3 
the  major  effort  in  this  investigation.   Because  of  the  signal  attenuation  of  about  10 

when  using  a  diffuser,  measurements  will  probably  not  be  practical  below  about  292  or  293 

[nm] .   Even  down  to  these  wavelengths,  the  spectral  slit  width  will  have  to  be  increased  to 

2 

1  [nm]   and  the  entrance  slit  area  to  about  1  [mm  ].   Because  the  ratio  of  hemispherical 

spectral  irradiance  to  solar  spectral  irradiance  changes  orders  of  magnitude  less  relative 
to  wavelength  than  the  hemispherical  quantity  itself  [14],  the  direct  measurement  of  the 
ratio  should  be  considered.   If  the  need  for  going  to  shorter  wavelengths  is  sufficiently 
great,  development  of  a  lens  or  optical  system  capable  of  hemispherical  viewing  and  cosine 
law  response  might  be  attempted.  Even  if  this  were  successful  and  the  signal/noise  problem 
down  to  290  [nm]  were  solved,  the  difficult  problem  of  accurately  correcting  for  scat- 
tered flux  would  remain. 
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Finally,  any  extensive  research  in  UV  solar  terrestrial  radiometry  should  consider  an 
instrument  in  which  the  spectral  irradiances  at  all  wavelengths  of  interest  can  be  measured 
simultaneously.   Changes  in  atmospheric  conditions,  often  in  a  relatively  short  time,  make 
this  feature  highly  desirable.   Fourier  transform  spectroradiometers,  multi-channel  spec- 
troradiometers,  and  spectrographs  are  obvious  examples.   However,  in  order  to  make  accurate 
measurements  below  about   300  [nm] ,   this  instrument  must  have  a  slit-scattering  function 
comparable  to  that  of  a  high  quality  double  monochromator,  and,  to  the  best  of  our  knowl- 
edge, this  has  not  yet  been  realized. 

SUMMARY  of  CHAPTER  1.    This  chapter  addresses  the  difficult  problem  of  measuring  the 
solar  terrestrial  spectral  irradiance  from  about   290  [nm]  to  340  [nm] .    It  includes  an 
account  of  the  experimental  design  and  measurements  that  were  carried  out  in  Gainesville, 
Florida,  in  June  1980.   The  major  reasons  for  these  measurements  were  to   (1)  apply  current 
state-of-the-art  radiometric  instruments  and  techniques  to  UV  solar  terrestrial  measurements 
and  present  the  results  in  this  Manual  as  a  first  applications  chapter  and   (2)  explore 
the  potential  of  using  solar  data  at  wavelengths  below  305  [nm]   for  determining  atmos- 
pheric ozone  thicknesses.   Item  1  is  the  subject  of  this  chapter.   Item  2  is  being  developed 
for  submission  to  Applied  Optics. 

The  chapter  begins  with  estimates  of  the  spectral  irradiances  of  the  sun  (not  includ- 
ing sky)  as  seen  from  the  surface  of  the  earth.   These  are  calculated  using  the  equations 
developed  in  Appendix  1A.   The  results  of  the  calculations  are  shown  in  figures  2  through 
4.   Figure  3  shows  that  with  a  spectral  slit  width  of  0.003  [nm]   considerable  spectral 
structure  appears.   Therefore  if  spectral  slit  widths  of  0.1  [nm]  to  1.0  [nm]   are  used, 
as  is  common  in  radiometry,  the  measurements  will  provide  only  a  weighted  average  spectral 
irradiance,  weighted  by  the  slit-scattering  function  of  the  spectroradiometer. 

The  measurement  configuration  that  we  used  in  this  problem  is  shown  in  figures  5 
through  7.   The  advantages  of  the  common  practice  in  irradiance  measurements  of  using  a 
diffuser  as  the  first  element  of  the  spectroradiometer  are  reviewed.   However,  a  signal 
reduction  factor  of  about   1000  makes  the  diffuser  impractical  to  use  below  about   305 
[nm]   where  measurements  for  ozone  determinations  are  to  be  concentrated.   It  is  demon- 
strated that  if  a  diffuser  is  not  used,  the  flux  responsivity  function  must  be  sufficiently 
uniform  over  the  angular  extent  of  both  the  sun  and  the  lamp  standard  used  in  calibrating 
the  spectroradiometer. 

The  general  measurement  equation  for  incoherent  radiation  (p.  23,  910-3)  [9], 


S(Xo)   = 


J  J   (R00*L0  +  R01'L1  +  R02*L2  +  H03.L3).coBe.d».dA:4A       (1.3) 

^XJ AAJ Aw 

is  applied  to  viewing  both  the  sun  and  the  standard  lamp  resulting  in 
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where  Re(Ao,T)   is  the  irradiance  responsivlty  factor  at  wavelength   A   and  tracker 
orientation  T;   and  where  z   is  the  slit-scattering  function  (p.  15,  910-4)  [6]. 

The  estimated  irradiances,  measurement  configuration,  and  measurement  equation  are 
used  to  arrive  at  the  specifications  of  the  monochromator  required  for  these  measurements 
so  that  their  uncertainty  is  less  than  10  percent.   These  include  a  wavelength  uncertainty 

of  0.01  [rim]   or  less,  spectral  slit  widths  varying  from  about  0.2  to  1  [nm] ,   very  low 

—8 
(<10  )   wings  on  the  slit-scattering  function  and  as  large  a  geometrical  throughput  as 

available.   A  block  diagram  of  the  complete  spectroradiometer  selected  is  shown  in  figure 

5,  and  detailed  specifications  for  the  instrument  are  also  presented. 

In  order  to  minimize  errors  in  the  radiometric  measurements,  the  flux  responsivity  of 
the  spectroradiometer  is  characterized  for  the  parameters  for  which  a  significant  dif- 
ference exists  between  the  spectral  irradiance  of  the  sun  and  that  of  the  lamp  used  as  a 
calibration  standard.   The  parameters  involved  are  magnitude  of  the  flux  and  spectral 
distribution. 

For  magnitude,  the  following  correction  to  the  output  signal  is  experimentally 
obtained: 

|r  =  C(S*)   =   1  +  6.927-10"8-S'  +  7.174-10_15-S'2  (1.15b) 

where  S   is  the  "corrected"  output  signal  and  S'   is  the  output  signal  actually  observed. 

The  rapid  fall-off  in  solar  terrestrial  spectral  irradiance  below  305  [nm]   and  the 
extensive  structure  require  a  maximum  wavelength  uncertainty  of  about  0.01  [nm]   if  wave- 
length accuracy  is  not  to  be  a  significant  factor.   It  is  found  that  this  0.01  [nm] 
uncertainty  is  realizable  with  this  best  available  portable  instrument  only  by  using  the 
spectral  structure  of  the  sun  itself  in  a  wavelength  self-calibration  while  taking  solar 
data.   In  order  to  obtain  the  calibration  wavelengths  in  this  observed  spectral  structure 
(by  evaluating  eq.  (1.13)  as  a  function  of  wavelength)  the  central  portion  of  the  slit- 
scattering  function  is  required.   It  is  obtained  from  spectral  scans  of  various  mercury 
and  laser  lines. 

—8 
An  accurate  determination  of  the  distant  wings  (to  values  of  10    or  less)  of  the 

slit-scattering  function  is  required  in  order  to  perform  a  deconvolution  of  eq.  (1.13)  to 
obtain  a  correction  for  scattering.   This  has  not  been  possible  because  of  insufficient 
power   (<<  10  [mw])   in  available  lasers.   However,  an  estimate  of  the  required  long  wave- 
length wing  has  been  obtained  using  the  strong   253.7  [nm]   mercury  line.   The  result  is 

—8 
2*10    at   53  [nm]   from  the  peak,  with  an  uncertainty  of  about  a  factor  of  two. 
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The  detailed  manner  in  which  the  solar  measurements  were  made  is  presented.   This 
included  an  accurate  optical  alignment  using  a  laser  (required  because  a  diffuser  was  not 
used  for  most  of  the  measurements) ,  frequent  spectral  irradiance  and  wavelength  calibra- 
tions, selection  of  integration  times  and  total  time  for  experiments,  visual  monitoring, 
temperature  control,  weather  conditions  and  computer  usage,  including  control.   Four  types 
of  data  were  obtained.   These  consisted  of   (1)  sharp  ozone  cut-off  data   (A  <  295  [nm]), 
(2)  data  for  the  new-type  ozone  determinations   (295  <  A  <  305) ,   (3)  data  for  Dobson  type 
ozone  determinations   (305  <  X  <  340),   and   (4)  data  required  for  performing  convolutions. 
Examples  of  these  four  types  of  data  are  shown  in  tables  B-l  through  B-7  in  Appendix  IB. 

For  reducing  the  data  and  estimating  the  uncertainty  of  the  measurements,  the  measure- 
ment equation  is  reduced  to 


/ 


A  +0.22 
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E,(A,t)-z(A    ,A,T)-dA 
A  o 
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A   -0.22  C(S').S'(X    .t,T).E?J(A    ,t) 

rrcTr? s  KiK'V    =  — : — r — •      (1-22) 
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where  f 

BjjCX    .t.Ti) 
P     =     -f-2 (1.23) 

g 

is  the  polarization  correction  required  when  the  calibration  data  S'    are  obtained  with 

g 
tracker  orientation  T   and  solar  data  S'   with  tracker  orientation  T.    E^   is  the 

s  A 

spectral  irradiance  of  the  standard  lamp  as  used  in  the  calibration  and  t   is  the  trans- 
mittance  of  any  filters  inserted  in  the  optical  path  for  solar  measurements  but  not  inserted 
for  the  lamp  calibrations.   The  manner  of  obtaining  each  independent  variable  in  eq. 
(1.22)  and  its  uncertainty  is  discussed.   These  uncertainties  are  summarized  in  table  2. 
Measurements  made  at  wavelengths  greater  than  305  [nm] ,   where  it  is  possible  (with  good 
signal-to-noise  ratio)  to  use  a  diffuser,  showed  a  difference  of  5.3  percent  between 
using  and  not  using  a  diffuser.   The  diffuser  data  are  believed  to  be  more  accurate,  and 
therefore  all  non-diffuser  data  will  be  corrected  (multiplied  by  0.947).    This  correction 
has  not  yet  been  made  for  the  data  shown  in  tables  B-l  to  B-7  in  Appendix  IB. 

Spectral  scattering  is  found  to  be  insignificant  down  to  290  [nm] .    This  requires 
the  distant  wings  of  the  slit-scattering  function  to  be  about   10    .    It  is  believed 
that  this  value  was  realized  in  our  non-diffuser  measurements  because  the  irradiation  of 
the  slit  of  the  monochromator,  in  this  case,  was  quite  coherent. 
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The  uncertainty  associated  with  each  dependent  and  independent  variable  in  eq.  (1.22) 
is  discussed  and  estimated.   The  overall  or  total  uncertainty  is  obtained  by  taking  the 
square  root  of  the  sum  of  the  squares  of  the  individual  uncertainties  in  table  2.   This 
total  uncertainty  varies  from   7.4  percent   at  the  short  wavelengths  to   5.7  percent   at 
305  [nm]   and  longer  wavelengths. 

Comparison  of  a  number  of  experimental  runs  shows  that  the  repeatability  is  not  as 
good  as  the  uncertainty  estimates  predict.  An  additional  uncertainty  of  2  percent  is 
added  to  the  totals  given  above  to  take  account  of  this. 

Means  that  can  be  taken  to  reduce  these  uncertainties  are  discussed.   An  improvement 
of  about  one  third  is  considered  feasible  if  the  mirror  solar  tracker  can  be  eliminated 
and  the  monochromator  itself  mounted  on  a  precision  solar  tracker. 

In  some  applications  in  which  these  measurements  will  be  used,  only  ratios  of  solar 
terrestrial  spectral  irradiance  at  wavelength  intervals  of  only   2  or  3  [nm]   are  required. 
These  are  estimated  to  have  much  lower  uncertainties  —  about   3  percent. 

Various  suggestions  are  made  for  further  research  in  UV  solar  terrestrial  spectral 
irradiance  measurements.   These  include   (1)  direct  confirmation  of  the  effect  of  coherence 
on  scattering,   (2)  extension  of  the  solar  measurements  to  include  global  (sun  plus  sky) 
measurements,  and   (3)  the  development  of  an  instrument  in  which  data  for  all  wavelengths 
of  interest  are  obtained  virtually  simultaneously. 
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Appendix    1A    (Part    III) 

Computation  of  Solar  Terrestrial   Spectral 
Irradiance  between   280   and    310    [nm] * 


■^he  purpose  of   this  Appendix  is   to  present   our   currently  recommended  procedure   for   com- 
puting  the   solar   terrestrial   spectral   irradiance   in  the     280   to   310    [nm]      region.      It 
should  be  emphasized,    however,    that   the  accuracy  of   this   recommended  procedure  has  not  yet 
been  critically  evaluated. 
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In  irradiance  calculations,  as  in  most  calculations  involving  radiometric  quantities, 
there  is  less  chance  of  error  if  one  begins  with  the  quantity  radiance.   The  solar 
spectral  radiance  at  the  earth's  surface,   L  ,   is  equal  to  the  product  of  the  solar  spec- 

A 

tral  radiance  outside  the  earth's  atmosphere  (extraterrestrial)   L^   and  the  spectral 

A 

transmittance  of  the  atmosphere  x(X)   along  the  path  associated  with  this  radiance 
(p.  38,  910-1)  [15];  that  is 

L,   =   Lf»T(A)   [W-m"2-sr"1.nm"1].  (A-l) 

A       A 

The  two  spectral  radiances  in  eq.  (A-l)  are  shown  in  figure  A-l  along  with  the  major 
constituents  of  the  atmosphere  —  ozone,  aerosol,  and  air  —  that  absorb  or  scatter  the 
solar  flux  in  the  region  between  280  and  340  [nm] .    The  spectral  transmittance  asso- 
ciated with  each  of  these  constituents  is,  by  the  familiar  exponential  law  of  transmit- 
tance, 

-q. (X) -seq .9 
x.(X)   =  e  X  (A-2) 

where  q.(X)   is  the  optical  depth1  of  the   i    constituent  for  a  vertical  path  (zenith 
angle  0=0).   Seq.0   is  the  ratio  of  the  length  of  the  actual  path  of  the  direct  solar 

flux  to  that  of  the  vertical  path.   The  optical  depth  q.   can  also  be  written  as 

-3-1  1  2 

n.a.W.   or  k.W.   where  n.  [cm  -atm  is  the  number  density,  o.    [cm  ]  is  the  atten- 

11111          l                                     l  _i_i 

uation  (or  absorption)  cross  section,   W.  [cm«atm]   is  the  thickness,  and  k.  [cm   «atm  ] 

is  the  attenuation  (or  absorption)  coefficient.2  The  overall  transmittance,   x(X),   of 

the  atmospheric  model .shown  in  figure  A-l  is 


^he  usual  symbol  for  the  optical  depth  in  the  literature  on  atmospheric  optics  is  t. 
However,  we  have  already  used  x   throughout  the  rest  of  this  Manual  for  transmittance. 
Accordingly,  we  have  adopted  here  the  symbol  q  for  optical  depth,  writing  q.(X) 
instead  of  x.(X).    Incidentally,  don't  confuse  this  with  the  unit  of  photon  flux,  the 
quantum   [q],   which  we  always  enclose  in  square  brackets,  as  with  all  unit-dimension 
symbols. 

2The  CIE-IEC  International  Lighting  Vocabulary  (3rd  ed.)  calls  this  the  linear  attenuation 
(or  absorption)  coefficient. 
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SUN 


OZONE      (q3) 


AIR         (q.,) 


AEROSOL       (q2) 


EARTH'S  SURFACE 


Figure  A-l.    Simple  model  of  earth's  atmosphere  showing  the  major  constituents  that 
absorb  or  scatter  solar  flux  in  the  UVB  and  the  solar  terrestrial 

spectral  radiance  L   and  solar  extraterrestrial  spectral  radiance  L  , 

A  A 

q.   is  the  vertical  optical  depth  of  the  ith  constituent. 
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-q.-'seq-.U      -q    *seq„«      -q„«seq_0 
x(A)      =      t1(A).t2(A).t3(A)      =      e  -e  -e  (A-3) 

-(q1'seq16  +  q«seq28   +  q.-seq  J)  (A-4) 

t(A)      =      e 

and   substituting   eq.    (A-4)    into   eq.    (A-l) ,    we  obtain 

-(q1«seq19   +  q2«seq28   +  q^seq.^) 
LA     =      V6  '  (A"5) 

where  q1  ,  q„ ,   and  q_   represent  the  vertical  optical  depth  for  air,  aerosol  and  ozone, 
respectively. 

To  obtain  the  solar  terrestrial  spectral  irradiance  E.   from  the  solar  terrestrial 

A 

spectral  radiance  L^   we  apply  the  equation  (p.  12,  910-2)  [16] 


E   =   /  L  •cos6'  -do)  (A-6) 


to  the  configuration  shown  in  figure  (A-2) .   Since   9'   does  not  exceed  about  0.25  [deg] , 


/V 


E   =    /  L,  -dco  (A-7) 


U) 


to  a  very  high  accuracy.   Substituting  the  expression  for  L   in  eq.  (A-5)  into  eq.  (A-7), 
we  obtain 


-(q1«seq].e  +  q  -seq^  +  q3«seq39) 
E   =   |  L  «e  «dcj.  (A-8) 


"    / 


Since  each  of  the  q.   and  seq.  8  vary  only  slightly  over  m 


i 


-(q  -seq  6  +  q  -seq  G  +  q^seq^)     .  _2 

E,   =   e  •   /  L  -dco  [W«m   «nm   J.  (A-9) 

A  -'A 


•  /*. 


The  solid  angle  u>  subtended  by  the  sun  at  the  earth's  surface  is  insignificantly  different 
from  that  subtended  by  the  sun  at  the  top  of  the  earth's  atmosphere.   Therefore 

-(q  -seq  9  +  q,-seq  8  +  q_-seq  9) 
E   =   e   X    L     12-  J    J   ,Ee  (A_10) 

A  A 


55 


SUN 


APPROX. 
0.5  [DEG] 


SURFACE  OF  IRRADIATION 
NEAR  EARTH'S  SURFACE 


Figure  A-2.  Geometry  involved  in  relating,  through  eq.  (A-6) , 
the  solar  terrestrial  spectral  radiance  and  solar 
terrestrial  spectral  irradiance. 
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where  E.   is  the  solar  extraterrestrial  spectral  irradiance  and  E,   is  the  solar  terres- 

A  A 

trial  spectral  irradiance. 

The  expression  that  we  have  used  for  the  vertical  optical  depth  for  air  q..   is  [17] 

1.0456'P   /300\4   0.1462-(300/X)2  ,t        • 

ql  =    1013    '{—)    '   e  (A-U) 

where  P  is  the  atmospheric  pressure  in  millibars  [mb]  at  the  point  for  which  E.   is 

A 

being  determined  and   X   is  the  wavelength  in  nanometers  [nm] . 

The  vertical  optical  depth  for  aerosol  q   for  the  limited  range  of  290  to  340 
[nm]   may  be  expressed  [3]  as 

q2  =  q20  -  s20-(x-3oo)  (a-12) 

where  q    and  S„_   are  constants.   The  constants  used  in  the  computation  of  E   in 
figure  2  through  figure  4  were  representative  values  obtained  in  our  1980  Gainesville 
measurements   (q„n  =  0.5  and   S„n  =  0.0033  [nm  I).1 

The  vertical  optical  depth  for  ozone  q_  may  be  written 

q3(y)   =   k3(y)-W3(y)  (A-13) 

where  k_(y)  [cm   «atm   ]   is  the  ozone  absorption  coefficient  and  VL(y)  [cm*atm]   is  the 
ozone  vertical  thickness  at  a  point  at  an  altitude  y  [km]   where  E   is  being  computed. 

A 

The  absorption  coefficient   k~   is  also  a  function  of  temperature.   At  wavelengths  between 
245  [nm]   and  311  [nm] ,   which  includes  the  wavelengths  with  which  this  investigation  is 
primarily  concerned,  this  functional  dependence  is  not  great  and  may  be  represented  by  [17] 

k3(X,t)   =   k3(A,to).[l  +  |[(t-to)/642]  +  [(t-to)/395]2}-e(A-300)/4°]       (A-14) 

where  t  and  t   are  in  degrees  Celsius  [°C]   and   X   is  in  nanometers  [nm] .    For  our 
calculations,  values  of   k~(X,t  )   were  obtained  from  measurements  by  Bass  [18]. 

To  take  account  of  the  variation  of  k,   in  the  earth's  atmosphere  due  to  the  varia- 
tion of  ozone  temperature  and  number  density,  the  following  procedure  was  followed.   A 
differential  element  of  ozone  optical  depth  dq_   at  an  altitude  y  where  the  temperature 
is   t  may  be  written  as 


^hose  interested  in  how  q„   may  vary  for  different  types  of  aerosols  and  in  analytical 
expressions  covering  a  much  larger  wavelength  range  are  referred  to  Green  and  Schippnick 
[17]. 
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dq3   =   k3(t)-dW3(y).  (A-15) 

The  total  ozone  optical  depth  is 


•i 


W3(y) 


q3      =  k3(t)-dW3(y),  (A-16) 


where     W_(y)      is   the  thickness  of   ozone  above  altitude     y.      In   eq.    (A-16)   we  are   inte- 
grating  from  the  top  of   the  atmosphere      (y  =  °°)      where     W_   =   0     to   a  point    in  the  atmos- 
phere where  the  altitude   is     y     and   ozone  thickness   above     y      is     W.(y).      If     W  is  the 
total   thickness  of   ozone       W        =  W   (0)        and     N(y)      is   the   fraction  of   this   thickness  above 


altitude     y, 


dW3(y)      =     W3Q.dN(y)      =     W^-^Vdy  (A-17) 


and   from  eq.    (A-16) 

.N(y) 


q3(y)      =  k3(t)-W3Q-dN  (A-18) 


I 


q3(y)  k3(t).W30.p.dy.  (A-19) 


Substituting   eq.    (A-14)    into   eq.    (A-19),    we  obtain 


q3 


(y)     ■  /'[l  +  {[(t-t„)/M2]   +   [(t-to)/395]2}.e(*-300)/40].k3(to, 


•«30-f  K 


q3 


oo 

(y)      -      -k3(to).W30.J"[l+{[(t-to)/642]+[(t-to)/395]2}.e"-3('0»/4°] 


fK 
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(A-20) 


(A-21) 


The  total  ozone  optical  depth   q,.(0)   which  we  designate  simply   q„   is 


q3(0) 


-W^o-jJ 


■dy 


+  /{[<t-t0>/642]  +    [(t-to)/395]2}.e^-30°)^.(f).dy] 


(A-22) 


q3   =   -k3(to).W30- 


_  00  _ 

-1  +  ^-W">.j  |[(t-to)/642]+[(t-to)/395]2}.(f).dy  I.    (A-23) 


For  ozone,   N(y)   may  be  written  [19] 


N(y)      = 


0.1307 


0.9709 


2.350  +ey/2'66  1+e(y-22.51)/4.920 


(A-24) 


where  the  altitude  y  is  in  kilometers  [km].   Differentiating  eq.  (A-24),  we  obtain 


dN        0.04914-e 
dy 


y/2.66 


0.1973-e 


(y-22.51)/4.92 


(2.350  +e^/2-66)      (l  +  e(y-22.51)/4.92)' 


(A-25) 


and  substituting  this  into  eq.  (A-23) 


,  t^   -v  t,   [-,    ,    n      (X-300)/40| 
q3   =   k3(to).W30.[n-C.e^         J 


(A-26) 


where 


C     =    /  {[tt-V7642]   +   [Ct-^0)7395j2} 


,,.37/2.66  o.1973.e^-22-51)/4-92\ 

— nz  1  *dy.         (A-27) 


(2.350  +e^/2'66)2  (1+e(y-22.51)/4.92): 


The  constant  C  may  be  numerically  evaluated  for  a  particular  t   using  tables  giving 

atmospheric  temperature  t  as  a  function  of  altitude  y   [20] .   These  tables  depend  on  the 

latitude  and  the  time  of  year  (month).   For  T   =  18  [°C]   and  a  latitude  of   30°N   in 

o 
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July,   C  =  -0.0714.    For  the  same  latitude  but  the  month  of  January,   C  changes  only 
slightly  to  -0.0745.    For  T  =  25  [°C]   and  a  latitude  of   30°N   in  July,   C  =  -0.0762. 

For  A  <  310  [nm] ,   the  absorption  coefficient  for  ozone  at  a  particular  temperature 
t   may  be  approximated  by  the  expression  [17] 

k3(to'X)   =  Q/(A-300)/Y  ?  (A"28> 

p  +  e 

Performing  a  non-linear  least  squares  fit  of  Bass'  data  [18]  from  288  to  305  [nm]   to  this 
expression  results  in 

a  =  11.1419 
6  =   0.04928 
y  =   7.1231 

where  k~(t  ,A)   has  the  units  of   [cm   -atm  ]. 

Expressions  used  for  seq.6  . in  eq.  (A-10)  are  [21] 

seq^e.y.)   =  fl  -  ^_1  J  (A_29j 

where 

r.  =  [1  +  (yi/R)]2,  (A-30) 

and  where  0   is  the  solar  zenith  angle,   R  [km]   is  the  earth's  radius   (6371  [km])   and 
where  y.  [km]   is  the  effective  altitude  for  the   ith  species.   We  used  [21] 

S 

y  (air)   =  5.69  [km] 
y„ (aerosol)   =   1  [km] 
y„(ozone)   =   23.5  [km]. 

For  small  6,   seq, (0.,y.)  =  secant  6..    The  difference  between  seq.   and  secant   is  not 
ni   l'-'i  l  l 

significant  except  for  very  large  zenith  angles.   Even  at   0  =  58  [deg] ,   the  difference 

for  the  case  of  ozone  is  only  0.9  percent  which  amounts  to  a  difference  in  computed  E 

A 

from  eq.  (A-10)  of  about   1  percent   at   304.6  [nm]   and   3  percent  at   294.5  [nm] .    For 
air,  the  difference  between  using   seq0   and   secant  0   at   0  =  58  [deg]   is  about   0.2 
percent  which  amounts  to  about   1/2  percent   in  E   in  this  same  wavelength  range. 


60 


The  remaining  quantity  required  to  compute  the  terrestrial  solar  spectral  irradiance 
E^   is  the  extraterrestrial  solar  spectral  irradiance  Ee.   We  needed  high  spectral  reso- 
lution  (AA  <<  0.2  [nm])   data  for  E   in  order  to  determine  the  effect  of  instrumental 

spectral  slit  widths  equal  to  or  much  larger  than   0.2  [nm] .    High  resolution  E8   data 

A 
could  not  be  found  for  the  spectral  region  of  interest.   However,  high  resolution 

(AA  =  0.003  [nm])   L^  data  were  available  for  the  solar  center  and  one  portion  of  the 

solar  limb  [22].   These  L^  data  were  used  to  calculate  E8   in  the  following  manner.   We 

start  with  the  relationship  between  spectral  irradiance  and  radiance  in  eq.  (A-6) .   As 

before  [eq.  (A-7)],  for  the  case  of  terrestrial  solar  spectral  irradiance,  this  becomes 


/ 


L    •dto, 

A 


(A-7) 


SUN 


APPROX.  <(> 


L°«t» 


L6X(0) 


Figure  A-3.   Configuration  for  computing  solar  extraterrestrial  spectral 

e  e 

irradiance  E   from  solar  extraterrestrial  spectral  radiance  L  . 

A  A 


The  solid  angle  dw   subtended  at  the  earth's  surface  by  the  annular  ring  of  area 
2-n-r-a-d$      [see  figure  (A-3)]  on  the  sun's  surface  is 


2TT'r-a»d<j)'  costfr 
R2 


2TT»a    »sin(|)'d(j) -costj 
R2 


(A-31) 


where   R  is  the  distance  from  the  sun  to  the  earth.   Substituting  this   dw   in  eq.  (A-7) 
we  obtain 
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-(f)2-} 


it/2 
E®     =      2ir«  (-|)  •  I         L® -sincj)-  co s<f>  -d<j>.  (A-32) 


If   we   let      y    =   costf) ,      dy   =   -sin<J)'d<j)      and 


<   =   2 


EI   =   2 


2  r° 

17 '  (t)  '  I  LI^-(-d^)>  (A~33) 

2  r1 

*•  (f)  •  I  L^-ydy,  (A-34) 


3 

where  L  (y)   for   y  =  1   and   0.23   are  available  from  Kohl,  et  al  [22],   Following  Kohl 

A 

et  al  [23],   L  (y)   is  assumed  to  be  a  straight  line  from  y  =  1   to  y  =  0.23   and  from 

A 

y  =  0.23   to   y  =  0,   and  L(0)   is  assumed  to  be  zero.   This  simple  functional  dependence 
of  L  (y)   is  estimated  [23]  to  be  accurate  to  a  few  percent.   With  these  assumptions 

L^(0.23) 
LA(y)   =  — cT~23 —  *y   f°r  y  =  °  t0  °'23  (A-35) 

and 


e        L*(0.23)  -0.23L*(1)     (l*(1)  -L*(0.23))y 

L,  (v)   =   7T^ + 


AXH/  0.77  0.77 


(A-36) 


for   y  =  0.23   to   1. 


Substituting  this  into  eq.  (A-34)  and  performing  the  integrations,  we  obtain,  using 
a  =  6.96-108  [m]   and  R  =  1.496-1011  [m] , 

Ef  =  136-10_6-[0.2050.L^(0.23)  +  0.2861-L®"(1)1  .  (A-37) 

A  A  A 

Direct  measurements  of  (spectrally  averaged)  extraterrestrial  solar  irradiance  are 

available  for  a  spectral  slit  width  of  1.0  [nm]   [24]  and  0.34  [nm]   [25].   These  spec- 

e  e 

trally  averaged  quantities  E   are  related  to  E   by  the  expression 

A  A 
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I 


E^-z-dA 

£     =  ^ (A-38) 

•dA 


AA 
where   z   is  the  slit-scattering  function.   Computations  of   E   using  eqs.  (A-37)  and 

A 

(A-38)  and  a  triangular  slit  function  (best  estimate  available)  were  compared  to  the  direct 
measurements.   Though  the  relative  spectral  values  agreed  within  a  few  percent,  the  values 
themselves  differed  by  factors  of   0.852   (1.0  [nm]  slit)   and   0.813   (0.34  [nm]  slit), 
the  computed  being  larger  than  the  observed  in  each  case.   We  decided  to  use  the  more 
recent   1.0  [nm]   data  of  Heath  so  we  normalized  eq.  (A-37)  by  the  factor  of   0.852. 
That  is, 

E*:  =  0.852.136-10"6-[0.2050-L<T(0.23)  +  0. 2861-L<T(1)]  .  (A-39) 


Tables  A-l  and  A-2  give  values  of   E   convoluted  (eq.  A-38)  with  a  triangular  slit- 

A 

scattering  function  having  a  spectral  slit  width  AA   of   0.18  [nm]   and   1.0  [nm] , 
respectively. 

Finally,  we  have  all  that  is  needed  to  compute  the  terrestrial  solar  spectral  irradi- 

ance   E    for   A  <  310  [nm]   which  are  the  wavelengths  of  major  interest  to  us.   Figure  3 

-1 
is  the  result  of  this  calculation  using   q   =0.5   and  S„   =  0.003  [nm   ]   and  Bass' 

k„(t  ,A).    For  a  spectral  slit  width  AA  >>  0.003  [nm] ,   such  as  in  figures  2  or  4,  one 

must  convolute  the  computed   E,   with  the  slit-scattering  function  of  the  instrument;  that 

A 

is, 

/   E  -z-dA 

Aa 


Ear  =  — ^ (a-4o) 

/   z-dA 

-/AA 


For  figures  2  and  4,  triangular  slit-scattering  functions  were  used  with  the  spectral 
widths  indicated  in  the  figures. 
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ADDENDUM  to  Chapter  4  in  NBS  Tech  Note  910-2: 

p.  28;  1st  line  below  eq.  (4.35): 
Insert  a  superscript  one  -  -  -  inverse- square  law1  -  -  -  referring  to 
the  following  footnote  at  the  bottom  of  the  page: 


'See  Explanatory  Note  following  end  of  Appendix  4  on  p.  102, 


p.  102,  following  end  of  Appendix  4,  add: 

Explanatory  Note  concerning  unit-dimensions  in  eq.  (4.36)  on  p .  28  as  an 
expression  of  the  inverse-square  law. 

In  this  Manual  we  have  insisted  on  the  need  for  recognizing  and  treating  solid 

angles  as  having  a  unit-dimension  of   [sr],   as  in  Tables  4-1,  4-2,  and  4-3,  on  pp.  3  to 

6,  inclusive.   However,  in  eq.  (4.36)  there  appears  to  be  an  exception  to  this  where  we 

take  the  quotient  of   I  [W«sr   ]   and  the  square  of  D  [m]   and  equate  that  quotient  to 

_2 
dE  =  dF   [W*m   ],   seemingly  ignoring  the  unit-dimension   [sr].   Actually,  there  is  no 

inconsistency  and  the   [sr]   is  not  ignored.   To  demonstrate  this,  we  need  to  distinguish 
between  the  always  orthogonal  radial  lengths,  with  dimension   [L  ]   or  unit-dimension 
[m  ],   and  transverse  lengths,  with  dimension   [L  ]   or  unit -dimension   [m  ] ,   as  dis- 
cussed in  examining  the  concepts  of  plane  angle  and  solid  angle  and  their  dimensionality 
in  Appendix  2  of  NBS  Tech  Note  910-1,  particularly  on  p.  69.   In  those  terms,  the  unit- 
dimension  of  a  plane  angle  is   [6]  =  [m  *m   ]   and  that  of  a  solid  angle  is 

2   -2 
[to]  =  [m   *m  ]   (the  quotient  of  the  dimensions  of  an  intercepted  or  enclosed  area 

subtended  by  the  solid  angle  on  the  surface  of  a  sphere  about  its  vertex  divided  by  the 
dimensions  of  the  square  of  the  radius  of  that  sphere).  In  eq.  (4.36),  the  unit  dimen- 
sion of  D,   the  distance  along  the  direction  of  propagation,  is  clearly   [D]  =  [m  ] 

not   [ml.    But  that  of  the  normal  irradiance  dE   or  fluence  rate  dF  .   the  flux  per 
t  n  t 

_2 
unit  area  perpendicular  to  the  direction  of  propagation,  is   [dE  ]  =  [dF  ]  =  [W*m  ] 

_2 
not   [W*m  ] .    The  unit  dimensions  of  radiant  intensity  I  are  then 

—1        —2   2  2 

[I]  =  [W«sr   ]  =  [W*m   *m   ].    Accordingly,  the  unit  dimensions  of  the  quotient   I/D 

-2        -2   2   -2        -2 
are   [I'D  ]  =  [W*m  *m  *m  ]  =  [W*m  ]  =  [dF  ]  =  [dE  ],   so  we  are  consistent  after 

all. 
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